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When two solids are squeezed together they will in general not make atomic contact everywhere 
within the nominal (or apparent) contact area. This fact has huge practical implications and must 
be considered in many technological applications. In this paper I briefly review basic theories of 
contact mechanics. I consider in detail a recently developed contact mechanics theory. I derive 
boundary conditions for the stress probability distribution function for elastic, elastoplastic and 
adhesive contact between solids and present numerical results illustrating some aspects of the theory. 
I analyze contact problems for very smooth polymer (PMMA) and Pyrex glass surfaces prepared 
by cooling liquids of glassy materials from above the glass transition temperature. I show that the 
surface roughness which results from the frozen capillary waves can have a large influence on the 
contact between the solids. The analysis suggest a new explanation for puzzling experimental results 
[L. Bureau, T. Baumberger and C. Caroli, arXiv:cond-mat/0510232 about the dependence of the 
frictional shear stress on the load for contact between a glassy polymer lens and flat substrates. I 
discuss the possibility of testing the theory using numerical methods, e.g., finite element calculations. 
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FIG. 1: Two solid blocks squeezed in contact with the force 
F. The area of real (atomic) contact A is usually an extremely 
small fraction of the nominal (or apparent) contact area Aq. 

Appendix A: Power spectrum of frozen capil- 
lary waves 

Appendix B: Stress distribution in adhesive 
contact between smooth spherical surfaces 

1 Introduction 

How large is the area of real contact when two solids 
are brought into contact (fig. ^1? This fundamental ques- 
tion has interested scientists since the pioneering work of 
Hertz published 1882[l|. The problem is of tremendous 
practical importance, since the area of real contact in- 
fluence a large number of physical properties such as the 
contact resistivity, heat transfer, adhesion, friction and 
the wear between solids in stationary or sliding contact. 
The contact area and the gap between two solids is also 
of crucial important for seals and for the optical prop- 
erties of composite systems, e.g., the optical interference 
between glass lenses. 

Hertz studied the frictionless contact between elastic 
solids with smooth surface profiles which could be ap- 
proximated as parabolic close to the contact area. This 
theory predict that the contact area A increases non- 



linearly with the squeezing force Fn as A - fI'^. The 
simplest model of a randomly rough but nominal flat 
surface consist of a regular array of spherical bumps (or 
cups) with equal radius of curvature R and equal height. 
If such a surface is squeezed against an elastic solid with 
a flat surface, the Hertz contact theory can be (approxi- 
mately) applied to each asperity contact region and one 
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expect that the area of real contact should scale as 
with the load or squeezing force F^. However, this is not 
in accordance with experiment for randomly rough sur- 
faces, which shows that, even when the contact is purely 
elastic (i.e., no plastic flow), the area of real contact is 
proportional to F^ as long as the contact area A is small 
compared to the nominal contact area ^q- This is, for ex- 
ample, the physical origin of Coulombs friction law which 
states that the friction force is proportional to the normal 
load0. 

In a pioneering study, ArchardQ showed that in a 
more realistic model of rough surfaces, where the rough- 
ness is described by a hierarchical model consisting of 
small spherical bumps on top of larger spherical bums 
and so on, the area of real contact is proportional to the 
load. This model explain the basic physics in a clear 
manner, but it cannot be used in practical calculations 
because real surfaces cannot be represented by the ide- 
alized surface roughness assumed by Archard. A some- 
what more useful model from the point of view of appli- 
cations, was presented by Greenwood and Williamson 
They described the rough surface as consisting of spher- 
ical bumps of equal radius of curvature R, but with a 
Gaussian distribution of heights. This model predict that 
the area of real contact is nearly proportional to the load. 
A more refined model based on the same picture was de- 
veloped by Bush et al0. They approximated the sum- 
mits with paraboloids to which they applied the Hertzian 
contact theory. The height distribution was described by 
a random process, and they found that at low squeezing 
force the area of real contact increases linearly with F^. 
Several other contact theories are reviewed in Ref . |2] . 

All the contact mechanical models described above as- 
sumes that the area of real contact A is very small com- 
pared to the nominal contact area a-nd made up from 
many small Hertzian- like asperity contact regions. This 
formulation has two drawbacks. First it neglects the in- 
teraction between the different contact regions. That is, 
if an asperity is squeezed against a flat hard surface it will 
deform not just locally, but the elastic deformation field 
will extend some distance away from the asperity and 
hence influence the deformation of other asperities 
Secondly, the contact mechanics theories can only be ap- 
plied as long as A « Aq. 

A fundamentally different approach to contact me- 
chanics has recently been developed by Perssonj^ lioj . 
It starts from the opposite limit of very large contact, 
and is in fact exact in the limit of complete contact (that 



is, the pressure distribution at the interface is exact). 
It also accounts (in an approximate way) for the elas- 
tic coupling between asperity contact regions. For small 
squeezing force it predict that the contact area A is pro- 
portional to the load F-^, while as Fn increases A ap- 
proach ^0 in a continuous manner. Thus, the theory is 
(approximately) valid for all squeezing forces. In addi- 
tion, the theory is very flexible and can be applied to 
more complex situations such as contact mechanics in- 
volving viscoelastic solids0,Q|, adhesionfio] and plastic 
yield 9] . It was originally developed in the context of 
rubber friction on rough surfaces (e.g., road surfaces) 0, 
where the (velocity dependent) viscoelastic deformations 
of the rubber on many different length scales are accu- 
rately described by the theory. The only input related 
to the surface roughness is the surface roughness power 
spectrum C{q), which can be easily obtained from the 
surface height profile /i(x) (see Ref. 0|)- The height 
profile can be measured routinely on all relevant length 
scales using different methods such as the Atomic Force 
Microscope (AFM), and optical methodsfl^. 

In this paper I present a brief review of contact me- 
chanics. I also present new results and a few comments 
related to a recently developed contact mechanics the- 
ory. I derive boundary conditions for the stress proba- 
bility distribution function for elastic, elastoplastic and 
adhesive contact between solids and present numerical 
results illustrating some aspects of the theory. I analyze 
contact problems for very smooth polymer (PMMA) and 
Pyrex glass surfaces prepared by cooling liquids of glassy 
materials from above the glass transition temperature. I 
show that the surface roughness which results from the 
frozen capillary waves can have a large influence on the 
contact between the solids. I also discuss the possibility 
of testing the theory using numerical methods, e.g., finite 
element calculations. 

2 Surface roughness 

The most important property of a rough surface is the 
surface roughness power spectrum which is the Fourier 
transform of the height-height correlation function|l^: 

Here z = h{x.) is the height of the surface at the point 
X — (x, y) above a flat reference plane chosen so that 
(/i(x)) = 0. The angular bracket (...) stands for ensem- 
ble averaging. In what follows we consider only surfaces 
for which the statistical properties are isotropic so that 
C(q) = C{q) only depend on the magnitude |q| of the 
wavevector q. 

Many surfaces of interest are approximately self afline 
fractal. A fractal surface has the property that if one 
magnify a section of it, it looks the same. For a self afline 
fractal surface, the statistical properties of the surface are 



unchanged if we make a scale change, which in general is 
different for the perpendicular direction as compared to 
the parallel (in plane) directions. The power spectrum of 
a self affine fractal surface 

where the Hurst exponent H is related to the fractal di- 
mension via Df — 3 — H . In reality, surfaces cannot 
be self affine fractal over all length scales. Thus, the 
largest wave vector possible is gi ~ 27r/a where a is some 
atomic distance or lattice constant, and the smallest pos- 
sible wave vector is ql « 27r/L, where L is the linear size 
of the surface. Surfaces prepared by fracture of brittle 
materials may be self afHne fractal in the whole length 
scale regime from a to L, but most surfaces are (approxi- 
mately) self affine fractal only in some finite wave vector 
regime, or not self afhne fractal anywhere as for, e.g., sur- 
faces prepared by slowly cooling a glassy material below 
the glass transition temperature, see below. 

Randomly rough surfaces with any given power spec- 
trum C{q) can be generated usingp^ 



/i(x)=^i?(q)e'^ 



[q-x+0(q)] 



(1) 



where (/>(q) are independent random variables, uniformly 
distributed in the interval [0, 27r[, and with B{q) = 



1/2 

• 



(27r/L)[C(q)]i/2, where L = A, 

Technologically important surfaces can be prepared by 
blasting small hard particles (e.g., sand) against a smooth 
substrate. Such surfaces are typically (approximately) 
self affine fractal in some wave vector regime qo < q < qi 
as illustrated in Fig. [21 The roll-off wave vector go de- 
creases with increasing time of blasting. Sand blasted 
surfaces, and surfaces prepared by crack propagation, 
have typical fractal dimension Di w 2.2, corresponding 
to the slope —2{H -|- 1) = —3.6 of the logC — logg rela- 
tion for q > qQ. We now give two example of the power 
spectrum of rough surfaces. 

2.1 Road surfaces 

Asphalt and concrete road pavements have nearly 
perfect self-afRne fractal power spectra, with a very 
well-defined roll-off wave vector go — 2n/Xo of order 
1000 m~^, corresponding to Ao « 1 cm, which reflects 
the largest stone particles used in the asphalt. This is 
illustrated in Fig. 13 for two different asphalt pavements. 
From the slope of the curves for g > go one can de- 
duce the fractal dimension Di « 2.2, which is typical for 
asphalt and concrete road surfaces. The height distribu- 
tions of the two asphalt surfaces are shown in Fig.^ Note 
that the rms roughness amplitude of surface 2 is nearly 
twice as high as for surface 1. Nevertheless, the tire- 
rubber frictionfol is slightly higher on the road surface 
1 because it has slightly larger power spectra for most 
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FIG. 2; Surface roughness power spectrum of a surface which 
is self affine fractal for qi > q > qo. The long-distance roll- 
off wave vector go and the short distance cut-off wave vector 
qi depend on the system under consideration. The slope of 
the logC — logg relation for q > qo determines the fractal 
exponent of the surface. The lateral size L of the surface (or 
of the studied surface region) determines the smallest possible 
wave vector qL = 2n/L. 
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FIG. 3: The surface roughness power spectra C{q) for two 
asphalt road surfaces. 



g- values in Fig. U] Thus there is in general no direct cor- 
relation between the rms-roughness amplitude and the 
rubber friction on road surfaces. 

2.2 Surfaces with frozen capillary waves 

Many technological applications require surfaces of 
solids to be extremely smooth. For example, window 
glass has to be so smooth that no (or negligible) diffu- 
sive scattering of the light occur (a glass surface with 
strong roughness on the length scale of the light wave- 
length will appear white and non-transparent because of 
diffusive light scattering). For glassy materials, e.g., sili- 
cate glasses, or glassy polymers, e.g., Plexiglas, extremely 
flat surfaces can be prepared by cooling the liquid from 
well above the glass transition temperature Tg, since in 
the liquid state the surface tension tends to eliminate (or 
reduce) short-wavelength roughness. Thus, float glass 
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FIG. 4: The height distribution Ph for two different asphalt 
road surfaces. 



(e.g., window glass) is produced by letting melted glass 
flow into a bath of molten tin in a continuous ribbon. 
The glass and the tin do not mix and the contact sur- 
face between these two materials is (nearly) perfectly flat. 
Cooling of the melt below the glass transition tempera- 
ture produces a solid glass with extremely flat surfaces. 
Similarly, spherical particles with very smooth surfaces 
can be prepared by cooling liquid drops of glassy mate- 
rials below the glass transition temperature. Sometimes 
glassy objects are "fire polished" by exposing them to a 
flame or an intense laser beam^^ which melt a thin sur- 
face layer of the material, which will flow and form a very 
smooth surface in order to minimize the surface free en- 
ergy. In this way, small-wavelength roughness is reduced 
while the overall (macroscopic) shape of the solid object 
is unchanged. 

However, surfaces of glassy materials prepared by cool- 
ing the liquid from a temperature above Tg cannot be 
perfectly (molecularly) smooth, but exhibit a fundamen- 
tal minimum surface roughness, with a maximum height 
fluctuation amplitude of typically 10 nm, which cannot 
be eliminated by any changes in the cooling procedure. 
The reason is that at the surface of a liquid, fluctua- 
tions of vertical displacement are caused by thermally 
excited capillary waves (ripplons) J^il^]. At the surface 
of very viscous supercooled liquids near the glass tran- 
sition temperature these fluctuations become very slow 
and are finally frozen in at the glass transition jl5l|. 

The roughness derived from the frozen capillary waves 
is unimportant in many practical applications, e.g., in 
most optical applications as is vividly evident for glass 
windows. However, in other applications they may be 
of profound importance. Here we will show that they 
are of crucial importance in contact mechanics. For elas- 
tically hard solids such as silica glass, already a rough- 
ness of order 1 nm is enough to eliminate the adhesion, 
and reduce the contact area between two such surfaces 
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FIG. 5: The surface roughness power spectrum C(q) for a 
polished steel surface with the root-mean-square (rms) rough- 
ness 0.1 /im and the fractal dimension Dt — 2.2, and for a 
PMMA surface prepared by cooling the melted polymer be- 
low the glass transition temperature. In this case the rms 
roughness is 0.6 nm. 



to just a small fraction of the nominal contact area even 
for high squeezing pressures 0. This is the case even for 
elastically much more compliant solids such as Plexiglas 
(PMMA), see Sec. 10.1. 

The surface roughness power spectrum due to capillary 
waves is of the form (see Appendix A) [lil IitI Ii8| 



1 



(2) 



where 7 is the surface tension, k the bending stiffness, 
and p the mass density of the glassy melt. In Fig. [S] 
we compare the the surface roughness power spectrum 
C [q] for a (typical) highly polished steel surface with the 
root-mean-square (rms) roughness 0.1 p,va and the frac- 
tal dimension Df = 2.2, with that of a PMMA surface 
prepared by cooling the melted polymer below the glass 
transition temperature. In this case the rms roughness 
is 0.6 nm. Note that in spite of the much larger rms 
roughness of the polished steel surface, for very large 
wavevectors the power spectrum of the PMMA surface 
is larger. This fact is of extreme importance for contact 
mechanics, since the large g-region usually dominates the 
contact mechanics. 

In Sec. 10 we will study contact mechanics problems 
for PMMA and for Pyrex glass. The smallest relevant 
wavevector q = ~ 27r/L, where L is the diameter of 
the nominal contact area ^o- In the applications below 
L « 100 ^m (for PMMA) and L « l^m (for Pyrex) 
giving >> pg. Thus, the gravity term in (2) can be 
neglected. The mean of the square of the surface height 
fluctuation is given by 



{h') = f d\ C{q) = 27r r 

J J qi^ 



dq qC{q) 



27r7 



In 



qi_ 



( ^ 

27r7 V qL 



(3) 



where qc = (7/^)^/^ is a cross-over wavevector. Eq. 
(3) gives the rms roughness w 0.5 nm for PMMA and 
« 0.15 nm for Pyrex, which typicaUy correspond to max- 
imum roughness amphtude of about 5 nm and 1.5 nm, 
respectively (see Appendix A). 

3 Contact mechanics theories 

We first briefly review the Hertz contact theory for 
elastic spheres with perfectly smooth surfaces. Most con- 
tact mechanics theories for randomly rough surfaces ap- 
proximate the surface asperities with spherical or ellipti- 
cal "bumps" to which they apply the Hertz contact the- 
ory. In most cases the elastic coupling between the asper- 
ity contact regions is neglected which is a good approx- 
imation only if the (average) distance between nearby 
contact regions is large enough. A necessary (but not 
sufflcient) condition for this is that the normal (squeez- 
ing) force must be so small that the area of real contact 
is very small compared to the nominal contact area. This 
is the basic strategy both for the Greenwood- Williamson 
theory and the (more accurate) theory of Bush et al, and 
we describe briefly both these theories. The theory of 
Persson start from the opposite limit of complete contact, 
and the theory is (approximately) valid for all squeezing 
forces. 




FIG. 6: Contact between two elastic spherical objects. The 
spheres are squeezed together with the force F and a circular 
contact area (radius ro) is formed. 



3.1 Hertz theory 

The Hertz theory describes the contact between two 
elastic spherical bodies (radius Ri and R2, respectively) 
with perfectly smooth surfaces. Assume that the spheres 



are squeezed together by the force F. The deformation 
field in the solids can be determined by minimizing the 
elastic deformation energy. The radius rg of the circular 
contact region (see Fig. El is given by 



f = 



R1R2 
R1+R2 



1/3 



3f (1 - 



1/3 



where 



E 



El 



(4) 



(5) 



where Ei and E2 are the elastic modulus of the solids and 
vi and U2 the corresponding Poisson ratios. The distance 
s the two solids approach each other (the penetration, see 
Fig. O is given by 



R1R2 



/ 3F(1 



AE 



2/3 



(6) 



For the special case of a sphere (radius R) in contact with 
a flat surface we get from (5) and (6) the area of contact 



■kRs 



and the squeezing force 



AE 



3(1-^.2) 



,3/2^1/2 



(7) 



(8) 



The pressure distribution in the contact area depends 
only on the distance r from the center of the circular 
contact area: 



(T(r) = o-Q 



1/2 



where ctq = F/ tttq is the average pressure. 

The Hertz contact theory has been generalized to ad- 
hesive contact. In this case the deformation field in the 
solids is determined by minimizing the elastic energy plus 
the interfacial binding energy. In the simplest case it is 
assumed that the wall-wall interaction has an infinites- 
imal extent so that the binding energy is of the form 
7rrQA7, where A7 = 71 -f- 72 — 712 is the change in the 
surface free energy (per unit area) as two flat surfaces 
of the solids makes contact. The resulting Johnson- 
Kendall-Roberts (JKR) theory^! is in good agreement 
with experiment [20j. For solids which are elastically very 
hard it is no longer possible to neglect the finite extent 
of the wall-wall interaction and other contact mechan- 
ics theories must be used in this case. The limiting case 
of a rigid sphere in contact with a flat rigid substrate 
was considered by DeriaguinEDl in 1934, and result in a 
pull-off force Fpuii-off — 27TA'yR which is a factor of 4/3 
larger as predicted by the JKR theory. A more accurate 
theory, which takes into account both the elasticity of the 



solids and the finite extent of the wall-wall interaction, 
will give a pull-off force somewhere between the JKR 
and Derjaguin limits. In Appendix B I discuss a partic- 
ular elegant and pedagogical approach to the adhesive 
contact mechanics for spherical bodies, which has been 
developed by Greenwood and Johnson, and which takes 
into account both the elasticity and the finite extent of 
the wall-wall interaction. 



FIG. 7: Contact between rough surface and a plane. All 
asperities with height z greater than the separation d make 
contact. 



and the nominal squeezing stress 




4S 



^0 3(1 - v^) 



no 



/ dh {h-df'^R^/^Ph. (12) 

Jd 



In Fig. |S1 1 show, for a typical case, the logarithm of 
the relative contact area A/Aq and the logarithm of the 
average area a of a contact patch, as a function of the log- 
arithm of the nominal squeezing pressure ao (in units of 
E/ [l — v"^)). Note that the contact area depends (weakly) 
non-linearly on the squeezing pressure (Tq and that the 
(average) contact patch area a increases with the squeez- 
ing pressure. Here we have defined 



where Ai is the area of contact patch i. Another defini- 
tion, which gives similar results, is a = A/N = Ai/N. 



3.2 Greenwood- Williamson theory 

We consider the friction less contact between elastic 
solids with randomly rough surfaces. If z = ft.i(x) and 
/i2(x) describes the surface height profiles, and Ei and E2 
are the Young's elastic modulus of the two solids, and vi 
and V2 the corresponding Poisson ratios, then the elastic 
contact problem is equivalent to the contact between a 
rigid solid with the roughness profile /i(x) = /ii(x) -|- 
/i2(x), in contact with an elastic solid with a flat surface 
and with the Young's modulus E and Poisson ratio v 
chosen so that (5) is obeyed. 

Let us assume that there is roughness on a single length 
scale. We approximate the asperities as spherical bumps 
with equal radius of curvature R (see Fig. 0) ^ind with a 
Gaussian height distribution 

where h* is the root-mean-square amplitude of the sum- 
mit height fluctuations, which is slightly smaller than the 
root-mean-square roughness amplitude /irms of the sur- 
face profile ft,(x). Assume that we can neglect the (elas- 
tic) interaction between the asperity contact regions. If 
the separation between the surfaces is denoted by d, an 
asperity with height h > d will make contact with the 
plane and the penetration s = h — d. Thus, using the 
Hertz contact theory with s — h — d the (normalized) 
area of real contact 

AA r°° 

— — = nnoR / dh (h- d)Ph (10) 

where Aq is the nominal contact area and no the number 
of asperities per unit area. The number of contacting 
asperities per unit area 

— = no / dh Ph (11) 
^0 Jd 
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FIG. 8: The logarithm of the relative contact area A/Ao and 
the logarithm of the average area a of a contact region, as a 
function of the logarithm of the nominal squeezing pressure ctq 
(in units of E/{1 — v'^)). The results are for the Greenwood- 
Williamson model assuming roughness on a single length scale 
Am = 27r/gm with = 600 m^^ and ?m/irms = 1- 

If, as assumed above, the roughness occur on a single 
length scale, say with wavevectors with magnitudes in a 
narrow region around i.e., C{q) ~ Co(5((? — qm), then 
we can estimate the concentration of asperities, no, the 
(average) radius of curvature R of the asperities, and the 
rms summit height fluctuation h* as follows. We expect 
the asperities to form a hexagonal-like (but somewhat 
disordered) distribution with the "lattice constant" Am = 
271'/ (?m so that 

no « 2/iXl^3) - « 0.029<z^ 

The height profile along some axis x in the surface plane 
will oscillate ( "quasi-periodically" ) with the wavelength 
of order w Am roughly as h{x) « hiCOs{qmx) where hi = 



\/2/i,.ms- Thus, 

i « h"{0) = q^hi = y2gf„/i,„,. 

Finally, we expect the rms of the fluctuation in the sum- 
mit height to be somewhat smaller than /iq, typically 
h* « 0.6/irms- Nayakj^ has shown how to exactly relate 
no, R and h* to qm and h rins ■ 

Fuller and Tabor have generalized the Greenwood- 
Williamson theory to the adhesive contact between ran- 
domly rough surfaces!^. They applied the JKR asperity 
contact model to each asperity contact area, instead of 
the Hertz model used in the original theory. The model 
presented above can also be easily generalized to elasto- 
plastic contact. 

3.3 The theory of Bush, Gibson and Thomas 

The Greenwood- Williamson theory assumes roughness 
on a single length scale, which result in an area of real 
contact which depends (slightly) non-linearly on the load 
even for very small load. Bush et al0 have developed 
a more general and accurate contact mechanics theory 
where roughness is assumed to occur on many different 
length scale. This result in an area of real contact which 
is proportional to the squeezing force for small squeezing 
force, i.e., A ~ Fn as long as A « ^o- Thus, the area 
of real contact A (for small load) is strictly proportional 
to the load only when roughness occur on many different 
length scales. In this case the stress distribution at the 
interface, and the average size a of an asperity contact 
region, are independent of the applied load. That is, as 
the load is increased (but A << Aq) existing contact 
areas grow and new contact regions form in such a way 
that the quantities referred to above remains unchanged. 

Bush et al considered the elastic contact between 
an isotropically rough surface with a plane by approx- 
imating the summits of a random process model by 
paraboloids with the same principle curvature and ap- 
plying the Hertzian solution for their deformation. Like 
in the Greenwood- Williamson theory, the elastic inter- 
action between the asperity contact regions is neglected. 
For small load the theory predict that the area of real 
contact is proportional to the load, A = aF^^, where 

a^n^^^i^j d\q^C{q)^ (13) 

where n = {27^)^/^. 

3.4 Persson theory 

Recently, a contact mechanics theory has been devel- 
oped which is valid not only when the area of real con- 
tact is small compared to the nominal contact area, but 
which is particularly accurate when the squeezing force 
is so high that nearly complete contact occurs within the 
nominal contact area|3,@,^3- other contact mechan- 
ics theories 0, 0, |3 were developed for the case where 



C=ioo 




FIG. 9: An elastic (e.g., rubber) block (dotted area) in adhe- 
sive contact with a hard rough substrate (dashed area). The 
substrate has roughness on many different length scales and 
the rubber makes partial contact with the substrate on all 
length scales. When a contact area is studied at low magnifi- 
cation {( = 1) it appears as if complete contact occurs in the 
macro asperity contact regions, but when the magnification 
is increased it is observed that in reality only partial contact 
occurs. 

the area of real contact is much smaller than the nomi- 
nal contact area. The theory developed in Ref. 9, 10] 
can also be applied when the adhesional interaction is 
included. 

Fig.|51shows the contact between two solids at increas- 
ing magnification ^. At low magnification (C = 1) it 
looks as if complete contact occurs between the solids at 
many macro asperity contact regions, but when the mag- 
nification is increased smaller length scale roughness is 
detected, and it is observed that only partial contact oc- 
curs at the asperities. In fact, if there would be no short 
distance cut-off, the true contact area would vanish. In 
reality, however, a short distance cut-off will always exist 
since the shortest possible length is an atomic distance. 
In many cases the local pressure at asperity contact re- 
gions at high magnification will become so high that the 
material yields plastically before reaching the atomic di- 
mension. In these cases the size of the real contact area 
will be determined mainly by the yield stress of the solid. 

The magnification C refers to some (arbitrary) chosen 
reference length scale. This could be, e.g., the lateral size 
L of the nominal contact area in which case C, = L/ X — 
q/ qL, where A is the shortest wavelength roughness which 
can be resolved at magnification C^. Sometimes, when the 
surface roughness power spectrum is of the form shown 
in Fig. 121 we will instead use the roll-off wavelength Aq = 
27r/(7o as the reference length so that Q — Xq/X = q/q^. 

Let us define the stress distribution at the magnifica- 
tion C, 

P(a,C) = (<5(a-a(x,C))) (14) 
Here a(x, Q is the stress at the interface calculated when 



only the surface roughness components with wave vectors 
q < (ql is included. The angular brackets (...) in (14) 
stands for ensemble average or, what is in most cases 
equivalent, average over the surface area 



^0 J A 



(15) 



where A is the area of contact. If the integral in (15) 
would be over the whole surface area Aq then P{(t, () 
would have a delta function [{Aq — A)/Ao]6{a) but this 
term is excluded with the definition we use. The area of 
real contact, projected on the xy-plane, can be obtained 
directly from the stress distribution since from (15) it 
follows that 



f)^/..P(.,C) 



(16) 



We win often denote A{C)/Ao = P(C). 

In Ref. I derived an expression for the stress dis- 
tribution P(cr, C) at the interface when the interface is 
studied at the magnification C = L/X, where L is the 
diameter of the nominal contact area between the solids 
and A the shortest surface roughness wavelength which 
can be detected at the resolution C^. That is, the surface 
roughness profile at the magnification ^ is obtained by re- 
stricting the sum in (1) to q < QlC- Assuming complete 
contact one can show that P{(t, C) satisfies the diffusion- 
like equation 



where 



/(C) = 



E 



(17) 



(18) 



where qL — 2tt/L and q = (q^. It is now assumed that 
(17) holds locally also when only partial contact occur at 
the interface. The calculation of the stress distribution 
in the latter case involves solving (17) with appropriate 
boundary conditions. 

If a rectangular elastic block is squeezed against the 
substrate with the (uniform) stress CTo i then at the lowest 
magnification C = 1 where the substrate appear flat, we 
have 



P{a, 1) = Sia - ao) 



(19) 



which form an "initial" condition. In addition, two 
boundary conditions along the cr-axis are necessary in or- 
der to solve (17). For elastic contact P{(J, C) must vanish 
as (7 — > oo. In the absence of an adhesional interaction, 
the stress distribution must also vanish for a < since 
no tensile stress is possible without adhesion. In fact, we 
will show below (Sec. 4) that P{a, C) must vanish con- 
tinuously as fj — !■ 0. We will also derive the boundary 



conditions for elastoplastic contact (Sec. 6 and 7), and 
for elastic contact with adhesion (Sec. 8). 

Eq. (17) is easy to solve with the "initial" condition 
(19) and the boundary condition P(0, C) = 0. The area 
of (apparent) contact when the system is studied at the 
magnification ^ is given by 



Ao 



— / dx e-" = erf (1/2VG) (20) 
Jo 



where 



G(C) 



E 



(1 - iy^)<yo 



C<?L 



dq q^Ciq) 



(21) 



In Fig. HJa) we show the dependence of the normalized 
contact area (at the highest magnification) A/Aq on the 
squeezing pressure ao (in units of E / {1 — v'^)). In (b) we 
show how A{Q/Ao depend on the magnification (where 
C = 1 now refers to the resolution Aq — 2n/qo rather than 
the size L of the surface) for cro(l - ^'^)/E = 0.001. The 
results have been obtained from (20) and (21) for the 
surface "polished steel" with the power spectra shown 
shown in Fig. El 

When the squeezing force Fn — (TqAo is so small that 
A << Ao, the equations above reduces to A = aFj^ with 
a given by (13) with k — (S/tt)^/^, i.e., the Persson the- 
ory also predict that the contact area increases linearly 
with the squeezing force Fn as long as the contact area 
is small compared to the nominal contact area. 

3.5 Comparison with numerical results 

The theory of Bush et al, and of Persson predict that 
the area of contact increases linearly with the load for 
small load. In the standard theory of Greenwood and 
Williamson 4] this result holds only approximately and 
a comparison of the prediction of their theory with the 
present theory is therefore difficult. The theory of Bush 
et al.|^ includes roughness on many different length 
scales, and at small load the area of contact depends 
linearly on the load A = q;-Fn, where 



(1 



a 



E 



d'q q'C{q) 



-1/2 



(22) 



where n = (27r)^/^. The same expression is valid in the 
Persson theory but with n — (B/tt)^/^. Thus the latter 
theory predict that the contact area is a factor of 2/7r 
smaller than the theory of Bush et al. Both the theory of 
Greenwood and Williamson and that of Bush et al. , as- 
sume that the asperity contact regions are independent. 
However, as discussed in Ref. 3, for real surfaces (which 
always have surface roughness on many different length 
scales) this will never be the case even at a very low 
nominal contact pressure. As we arguedj^ this may be 
the origin of the 2 / vr-difference between the Persson the- 
ory (which assumes roughness on many different length 
scales) and the result of Bush et al. 
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FIG. 10: (a) The dependence of the normahzed contact area 
(at the highest magnification) A/Aq on the squeezing pressure 
(Jo (in units oi E/{1 ~ u^)), and (b) on the magnification for 
(7o(l - v^)IE = 0.001. For the surface "polished steel" with 
the power spectra shown shown in Fig. |3 



The predictions of the contact theories of Bush et al.p 
and Perssonlsl were compared to numerical calculations 
(see Ref. j^jlal). Borri-Brunetto et al. sT] studied the 
contact between self affine fractal surfaces using an es- 
sentially exact numerical method. They found that the 
contact area is proportional to the squeezing force for 
small squeezing forces. Furthermore, it was found that 
the slope a(C) of the line A ~ a(C)-FN decreased with 
increasing magnification ^. This is also predicted by the 
analytical theory [Eq. (22)]. In fact, good agreement 
was found between theory and computer simulations for 
the change in the slope with magnification and its depen- 
dence on the fractal dimension ZJf . 

Hyun et al. performed a finite-element analysis of con- 
tact between elastic self-afhne surfaces. The simulations 
were done for a rough elastic surface contacting a per- 
fectly rigid flat surface. The elastic solid was discretized 
into blocks and the surface nodes form a square grid. The 
contact algorithm identified all nodes on the top surface 
that attempt to penetrate the flat bottom surface. The 
total contact area A was obtained by multiplying the 
number of penetrating nodes by the area of each square 
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FIG. 11: Dots: the factor function of Hurst's exponent 

H for self affine fractal surfaces. The two horizontal lines are 
the predictions of the theories of Bush et al. (solid line) and 
Persson (dashed line). From Ref. |25|. 

associated with each node. As long as the squeezing force 
was so small that the contact area remained below 10% 
of the nominal contact area, A/Aq < 0.1, the area of real 
contact was found to be proportional to the squeezing 
force. Fig. 1111 shows Hyun et al.'s results for the factor 
K in (22) as a function of Hurst's exponent H for self 
affine fractal surfaces. The two horizontal lines are pre- 
dictions of the theories of Bush et al. (solid line) and 
Persson (dashed line). The agreement with the analyti- 
cal predictions is quite good considering the ambiguities 
in the discretization of the surface (Sec. 12). The al- 
gorithm only consider nodal heights and assumes that 
contact of a node implies contact over the entire corre- 
sponding square. This procedure would be accurate if 
the spacing between nodes where much smaller than the 
typical size of asperity contacts. However, the majority 
of the contact area consists of clusters containing only 
one or a few nodes. Since the number of large clusters 
grows as H ^ 1, this may explain why the numerical 
results approach Persson's prediction in this limit. 

Hyun et al. also studied the distribution of connected 
contact regions and the contact morphology. In addition, 
the interfacial stress distribution was considered, and it 
was found that the stress distribution remained non-zero 
as the stress cr — > 0. This violates the boundary condition 
that P(cr, C) = for (T = 0. However, it has been shown 
analytically^] that for "smooth" surface roughness this 
latter condition must be satisfied, and the violation of 
this boundary condition in the numerical simulations re- 
flect the way the solid was discretized and the way the 
contact area was defined in Hyun et al.'s numerical pro- 
cedure (see also Sec. 12). 

In Ref. |53| the constant a was compared to (atom- 



istic) Molecular Dynamics (MD) calculations for self 
affine fractal surfaces (with fractal dimension Df = 2.2). 
The MD-result for a was found to be « 14% larger 
than the theoretical prediction, which is very similar to 
the numerical results of Hyun et al^^. Thus, for the 
same fractal dimension Di — 2.2 as in the MD calcu- 
lation, Hyun found « 13% larger a than predicted by 
the analytical theory (see Fig. Ill|l . Very recently, M. 
Borri-Brunetto (private communication) has performed 
another discretized continuum mechanics (DCM) calcu- 
lation (for exactly the same surface roughness profile as 
used in the MD calculation) using the approach described 
in Ref. and found a (at the highest magnification) 
to be « 20% smaller than the theoretical value. Tak- 
ing into account the atomistic nature of the MD model 
(and the resulting problem with how to define the con- 
tact area), and the problem with the finite grid size in 
the DCM calculations, the agreement with the theory of 
Persson is remarkable good. Since the Persson theory 
becomes exact as one approach complete contact, it is 
likely that the theory is highly accurate for all squeezing 
forces Fn- 




FIG. 13: The logarithm of the (normalized) contact area as 
a function of the logarithm of the magnification. The contact 
regions at the magnification (m , where the second derivative of 
the curve is maximal, is denoted as the macroasperity contact 
regions. The elastic modulus E — 10^^ Pa and Poisson ratio 
v — 0.3. For the surface "pohshed steel" with the power 
spectra shown in Fig. |3 and with the (nominal) squeezing 
pressure ctq = 10 MPa. 
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FIG. 12: The contact area at increasing magnification (a)- 
(d). The macro-asperity contact area (c) breaks up into 
smaller contact areas (d) as the magnification is increased. 



3.6 Macroasperity contact 

Consider the contact between two solids with nominal 
flat surfaces. Assume that the surface roughness power 
spectrum has the qualitative form shown in Fig. |21with 
a roll-off wavevector qq corresponding to the magnifica- 
tion C = 1. Fig. 1121 shows schematically the contact 
region at increasing magnification (a)-(d). At low mag- 
nification C < 1 it appears as if complete contact occur 
between the solids. At the magnification C ~ 1 a few 
non-contact islands will be observed. As the magnifica- 
tion is increased further the non-contact area will perco- 
late giving rise to contact islands completely surrounded 



by non-contact area. 

Fig. ^] shows, for a typical case, the logarithm of 
the (normalized) contact area as a function of the log- 
arithm of the magnification. We denote the contact re- 
gions at the magnification Cm, where the second deriva- 
tive of the curve is maximal, as the macroasperity con- 
tact regions. The macroasperity contact regions typi- 
cally appear for a magnification somewhere in the range 
Cm = Qm/Qo ~ 1 — 5, where A/Aq is below the (con- 
tact) percolation threshold (which typically occur when 
A/Aq « 0.3). The macroasperity contact regions [fig. 
I12r c)] breaks up into smaller contact regions (d) as the 
magnification is increased. 

The concept of macroasperity regions is very important 
in rubber friction, in particular in the description of the 
influence of the flash temperature on rubber friction[23j. 

3.7 On the validity of contact mechanics theo- 
ries 

The Greenwood- Williamson theory and the theory of 
Bush et al (and most other contact theories) are based 
on the assumption that the (average) separation between 
nearby contact regions is large enough so that the elastic 
deformation field arising from one contact region has a 
negligible influence on the elastic deformation at all other 
contact regions. This assumption is strictly valid only at 
low enough squeezing forces and for the magnifications 
C ~ Cm, where the contact regions takes the form shown 
schematically in Fig. I12r c). The reason for why, in gen- 
eral, one cannot apply these theories at higher magnifica- 
tion is that for magnification C > Cm, the macro-asperity 
contact regions break up into smaller closely spaced as- 
perity contact regions (see Fig. I12|l and it is, in gen- 



eral, not possible to neglect the elastic coupling between 
these microasperity contact regions. We note that the 
(average) pressure in the macroasperity contact regions 
is independent of the squeezing pressure ctq (but the total 
macroasperity contact area is proportional to ctq) so that 
the failure of the Greenwood- Williamson and Bush et al 
theories to describe the contact mechanics at high mag- 
nification is independent of the nominal contact pressure. 
The theory of Persson does (in an approximate way) ac- 
count for the elastic coupling, and may hold (approxi- 
mately) for all magnifications. 

4 Elastic contact mechanics 

Consider an elastic solid with a nominal flat surfaces 
in contact with a rigid, randomly rough substrate. We 
assume that there is no adhesive interaction between the 
solids. In this case no tensile stress can occur anywhere 
on the surfaces of the solids, and the stress distribution 
P(cr, C) must vanish for cr < 0. For this case one can 
show that P{a; vanish as cr ^ 0. This can be proved 
in two ways. First, it is known from continuum mechanics 
that when two solids with smooth curved surfaces are in 
contact (without adhesion) the stress in a contact area 
varies as ^^^^ with the distance ^ from the boundary line 
of the contact area. From this follows that P{o,C) ~ cr as 
cr — > and in particular P(0, C) = 0, see Ref. ^ . For the 
special case of bodies with quadratic surface profiles the 
stress distribution is known analytically. In particular, 
for spherical surfaces the contact area is circular (radius 
ro) and the pressure distribution 



a{r) = cr„ 



1/2 



where Cmax = 3cro/2, where cro is the average normal 
stress in the contact area. Thus, in this case 



1 



d^x S{a — cr(r)) 



2cr 



(23) 



which vanishes linearly with cr. 

A second proof follows from (17), and the observation 
that at any magnification the contact pressure integrated 
over the (apparent) area of contact must be equal to the 
load. Thus, multiplying (17) with a and integrating over 
a gives 



d 



oo Q2p 



-/(C) 



da— (a,C)=/(C)P(0,C) 



(24) 



where we have performed a partial integration and as- 
sumed that adP{a,C)/d<T ^ as cr —> oo (which holds 



because P(cr, () decreases exponentially or faster with in- 
creasing cr, for large cr). Since 

/■oo 

/ da aP{<7, C) ctq 
JQ 

for all magnifications we get 



del ^--^(-'0 = 



(25) 



Thus, since /(C) ^ we get from (24) and (25) that 
P(0, C) — 0. Hence, (24) together with the observation 
that the load must be carried by the area of (apparent) 
contact at any magnification, imply P(0, C) — 0. The fact 
that this (exact) result follows from (24) is very satisfac- 
tory and is not trivial since equation (24) is not exact 
when partial contact occur at the interface. Thus, the 
fact that (24) is consistent with the correct boundary 
condition lend support for the usefulness of the present 
contact mechanics theory. 
Integrating (24) over cr gives 



/>oo />C 

/ daP{<j,C) = / 

jQ Jo 



d_ 

and since 



da /(C) 



d^P 
5^ 



-/(C)|^(o,C) 



r da P{aX)^^^P{a) 
Jo ^0 



we get 



A' = ~Ao/(C)|^(0,C) 



where the index "el" indicate that the solids deform elas- 
tically in the area of contact. Since the total area (pro- 
jected on the xy-plane) is conserved A^i + Anon — Aq, 
where Anon denote the non-contact area, we get A'^^ = 



-A'„„ so that 



<o„ = ^o/(C)|^(o,C) 



(26) 



In analogy with heat current, we can interpret this equa- 
tion as describing the flow of area from contact into non- 
contact with increasing magnification. 

5 Elastic contact mechanics with graded elas- 
ticity 

In many applications the elastic properties of solids 
are not constant in space but varies with the distance z 
from the surface into the solid. Thus, for example, many 
solids have modified surface layers, e.g., thin oxide layers 
with different elastic properties than the bulk. In these 
cases the elastic modulus E = E{z) will depend on the 
distance z from the surface. The same is true for many 
biological systems where natural selection has optimized 
the surface properties for maximum functionality. Thus, 




FIG. 14: Vertical cut through the adhesion pad ol cicada. 
The foam-like structure of the adhesion pad makes it elasti- 
cally soft, which is necessary for the pad to adhere to rough 
surfaces. In addition, the elastic properties are graded |2a| 
which may also be important for the adhesion (see Sec. 8.5). 
Adapted from Ref. i27| . 



6 Elastoplastic contact mechanics: constant 
hardness 

We assume ideal elastoplastic solids. That is, the solids 
deform elastically until the normal stress reach the yield 
stress cry (or rather the indentation hardness, which is 
~ 3 times the yield stress of the material in simple com- 
pression), at which point plastic flow occur. The stress 
in the plastically deformed contact area is assumed to 
be (Ty (e.g., no strain hardening). Contact mechanics 
for elastoplastic solids was studied in Ref. "o), where 
the stress distribution was calculated from (17) using the 
boundary conditions P(0, C) = and P{<ty, C) = 0. The 
first (elastic detachment) boundary condition was proved 
above and here we show that the second (plastic yield) 
boundary condition follows from Eq. (17) together with 
the condition that at any magnification the contact pres- 
sure integrated over the (apparent) area of contact must 
be equal to the load. In this case the load is carried in 
part by the plastically deformed contact area. 

Integrating (17) from ct = to cry gives 



for example, the adhesion pads (see Fig. I14|l on many 
insects consist of graded materials which may exhibit ro- 
bust fiaw tolerant adhesion (see Sec. 8.5). It is therefore 
of great interest to study contact mechanics, both with 
and without adhesion, for solids with elastic properties 
which depend on the distance from the surface into the 
solid. 

The theory developed in Ref.0, 0| can be (approxi- 
mately) applied also to study contact mechanics for solids 
with graded elastic properties. Thus, the deformation 
field in a solid in a response to a surface stress distribu- 
tion which varies parallel to the surface with the wave- 
length A (wave vector q = 2tt/X) extend into the solid 
a distance of order ~ and is approximately of the 
form ~ exp(— gz). Thus, if we replace E in (18) with 
E{C) defined by 

f„°° dz E(z)exp(~qz) r°° .... 
E{C) = -^\oo / ^ , , ' =q dz E{z)cxp{~qz) 
Jp dz exp(-gz) Jo 

where q = qo(, we account for the fact that the effec- 
tive elastic modulus depends on the wavelength of the 
roughness profile. This approach should work best if the 
spatial variation in E{z) is weak, e.g., E'{z)/E{z) < q 
for z ~ 1/q. The same approach can be used when adhe- 
sion is included, but then the elastic modulus which enter 
in the detachment stress, and in the expression for the 
elastic energy, must also be replaced by the expression 
E{C) defined above. The limiting case of a thin elastic 
plate can be treated using the formalism developed in 
Ref. SEi. 



= -/(C)§^(0,C) + /(C)|^(ay,C) (27) 

Since 

'"da P(a,C)=P(C)-^ 

^0 

where Ad is the "elastic contact area" , we can write (27) 
as 

Ai,(C) = -Ao/(C)|^(0, C) + Aof{0^{<TYX) 



= -<o„(C) - 4i(C) 

where we have used (26), and where 



da 



(28) 



(29) 



is the contact area where plastic deformation has oc- 
curred. Thus, (28) is just a statement about conservation 
of area (projected on the xy-plane): 

^ci(C)+^;i(C) + <o„(C) = o, 

or after integration 

Ael + Api + Anon = Aq 

Next, let us multiply (17) with a and integrate over tr: 

8-cl ^^<^noX)-l rf-/(cy^ 



= /(CW|^(aY, 0-/(0 



dP 

da^(-,0 



/(O 



(30) 



Since 



^^ei(0 =^0 / "d^T <yP{<jX) 

Jo 

is the normal force carried by the elastically deformed 
contact area, (30) gives 

dP 



Using (29) this gives 

^^,'i(0 = -^YA;,i-Ao/(OP(aY,0 



(31) 



Since the normal stress in the contact region where plas- 
tic yield has occurred is equal to cry it follows that the 
plastically yielded contact regions will carry the load 
^pi — CTy^pi- Thus, since cty is assumed independent 
of C we get i^pj = fY^pi and (31) can be written as 



^ci(0 = -^;i(0-^o/(0^'('TY,0 



(32) 



But the total load F-^ = ctq^o — Fc\{C) + -Ppi(0 niust be 
independent of the magnification so that F^^{()+F^^{C) = 
0. Comparing this with (32) gives the boundary condi- 
tion P(cty,0 = 0. 

In Ref. |g| the "diffusion equation" (17) was solved 
analytically with the boundary conditions P{Q, C) = 
and P{aY,C) = 0- Qualitative we may state that with 
increasing magnification the contact area diffuses over 
the a — boundary into no-contact, and over the a = ay 
boundary into plastic contact. 

As an illustration Fig. 1 151 shows how (a) the logarithm 
of the (normalized) elastic contact area, and (b) the (nor- 
malized) plastic contact area depend on the logarithm 
of the magnification. The calculation is for the surface 
"polished steel" with the power spectrum shown in Fig. [3 
and with the (nominal) squeezing pressure uo = 10 MPa. 
The elastic modulus E — 10^^ Pa, Poisson ratio ly = 0.3 
and for several yield stresses indicated in the figure. Note 
that with increasing magnification, the plastically yielded 
contact area Api/Aa increases continuously towards the 
limiting value CTo/cry. 

7 Elastoplastic contact mechanics: size- 
dependent hardness 

Many materials exhibit a penetration hardness which 
depend on the size of the indentation. That is, cry is a 
function of In this case 
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FIG. 15: The dependence of (a) the logarithm of the (nor- 
malized) elastic contact area and (b) the (normalized) plastic 
contact area on the logarithm of the magnification. The elas- 
tic modulus E = 10^^ Pa, Poisson ratio v — 0.3 and for 
several yield stresses indicated in the figure. For the surface 
"polished steel" with the power spectra shown in Fig. |^and 
with the (nominal) squeezing pressure ao = 10 MPa. 



= a'^{OP{aYX) - /(0|^(0, + fiO^i^vX) 

which differ from (27) by the term crY(C)^'(o'y , C)- In this 
case area conservation require the definition 

A'^.iO = -A,ai,{OP{aY,0~Aof{0^{aY,0 (33) 

A similar "extra" term will appear in (30) but (31) will 
remain unchanged. Thus, since i^pi(C) = o'y(0^pi(C) + 
^y(0^pi(0, (31) takes the form 

^ci(0 = -^;i(C) + <tV(0^pi(0 - Ao/(OP(ay, 
Hence, in this case 

aV(OApi(0 - Ao/(OP(ay, = (34) 

This equation together (33) gives the correct boundary 
condition when cry(C) depend on 

8 Elastic contact mechanics with adhesion 

In this section I first derive and discuss the boundary 
conditions for the stress distribution relevant for adhesive 



contact between solids. I also derive an improved expres- 
sion for the elastic energy stored at the interface in an 
adhesive contact. Numerical results are presented to il- 
lustrate the theory. Finally, I present a few comments 
related to flaw tolerant adhesion, adhesion involving vis- 
coelastic solids, and adhesion in biological systems. 

8.1 Boundary condition and physical interpre- 
tation 

When adhesion is included, the stress distribution 
function P{<7, C) must vanish for a < —a^, where cra(C) is 
the highest tensile stress possible at the interface when 
the system is studied at the magnification We will now 
prove the boundary condition P{—(Jii{(), () = 0. When 
the adhesion is neglected da = and P(0, C) = which 
is the limiting case studied in Sec. 4. 

In the present case, the (normalized) area of (apparent) 
contact when the system is studied at the resolution is 
given by 



^(0 = 



MO 

An 



daP{a,0 



(35) 



where Aq is the nominal contact area. The area of real 
contact is obtained from (35) at the highest (atomic) res- 
olution, corresponding to the magnification Qi = L/a, 
where a is an atomic distance. 
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theory of cracks which gives (see Fig. I16|l 



V 1 - / 



(36) 



where a is a number of order unity (we use a — 9/8) 
and 7eff(C) is the (effective) interfacial binding energy 
when the system is studied at the magnification C,. At 
the highest magnification (corresponding to atomic reso- 
lution) 7eff = A7, where A7 = 71-1-72—712 is the change 
in the energy (per unit area) when two flat surfaces of the 
solids are brought into adhesive contact. The effective in- 
terfacial energy at the magnification C, is determined by 
the equation 



^*(C)7off(C)=^*(Ci)A7-t/ci(C) 



(37) 



where A*{C,) denote the total contact area which may 
differ from the projected contact area A{C,). C/ci(C) is the 
elastic energy stored at the interface due to the elastic 
deformation of the solid on length scales shorter than A = 
L/C, necessary in order to bring the solids into adhesive 
contact. 

Using (17) we get 



d_ 

dC J- 



— / da P{a, C) = a;(C)P(-aa, C)+ / da /(C) 



d^P 
9^ 



= CTl(C)P(-aa,C)-/(C)|^(-^a,C) 



or 



BP 

A'iO = Aoa',{C)P{-a,, C) - Aof{C)^{-<J,, C) (38) 
Next, using (17) we get 
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/oo 
da aP{aX) = -<(C)^a(C)^'(-'Ta, C) 




FIG. 16: An elastic solid (dotted area) in adhesive contact 
with a hard solid (black area) with surface roughness. The 
elastic solid will deform and fill out the substrate "cavity" only 
if the gain in adhesive energy ~ 7^^ is of order (or larger than) 
the elastic energy stored in the deformation field. The defor- 
mation field extend a distance d into the solid, and the strain 
is of order h/d giving the elastic energy of order d"^ E{h/d)^ . 
The detachment stress (Ja — Eh/d is defined as the stress 
when [/ad = Uei or yd « Eh^ so that Uad ~ {lE/dY^^. Since 
d is of order the wavelength of the surface roughness we get 
(36) to within a factor of order unity. 
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The detachment stress a^ can be determined using the 
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fTa7^(-fTa,C)+-P(-<^a,C) 

oa 



Thus, the normal force acting in the area of contact sat- 
isfies 

F'(C)--Ao<(CK(C)P(-aa,C) 



the stress distribution in the area of real contact. I find 
that the stress distribution is Hertzian-like but with a 
tail of the form ^ (ct + (Ta) extending to negative (ten- 
sile) stress. 



+^o/(C) 



Using (38) this gives 

F'iO = -^a(C)^'(C) + AoPi^a,, (39) 

Now, note that since A'{() < 0, both terms on the right 
hand side of (39) are positive. This implies that the nor- 
mal force F(Q acting in the area of contact cannot, in 
general, be equal to the external load Fjq = Aoao. This 
remarkable result is possible only if there is an attractive 
interaction between the non-contacting surfaces. How- 
ever, this is exactly what one expect in the present case 
for the following reason: The (effective) surface energy 
7eff is determined by the work to separate the surfaces. 
Since the stress field in our model is everywhere finite, 
in particular the detachment stress Ca is finite, the sur- 
face forces must extend over a finite extent in order for 
the work of separation to remain non-zero. Only when 
the stresses are allowed to become infinite is it possible 
to have a finite work of adhesion with contact interac- 
tion (i.e., interaction of infinitesimal extent) 7 = x 00. 
Since the detachment stress is finite in our model, the 
attractive interaction between the walls must have a fi- 
nite extent and exist also outside the area of contact, see 
Fig. El Thus, (39) cannot be used to determine the 
boundary condition at cr = —fa, but describes how the 
force acting in the contact area increases with increasing 
magnification, as a result of the adhesional "load" . How- 
ever, the boundary condition can be determined using 
the first method described in Sec. 4, by considering the 
stress distribution close to the boundary line between the 
contact and non-contact area, which we refer to as the 
detachment line. The detachment line can be considered 
as a crack edge, and from the general theory of cracks one 
expect the tensile stress to diverge as ^"^/^ with the dis- 
tance ^ from the detachment line. However, this is only 
true when the interaction potential between the surfaces 
is infinitesimally short ranged, which is not the case in 
the present study (see above), where the stress close to 
the boundary line (in the crack tip process zone) will be 
of the form ct = da + C^^/^, where ^ is the distance from 
the boundary line. This gives da/d£, ~ C"^/^ and fol- 
lowing the argument presented in Sec. 4 (and in Ref. 
[1|) this imply that the stress distribution will vanish lin- 
early with (7 -I- (Ta as (T ^ — (Ta so that P( — (Ta, C) = 0. 

For the special case of a circular asperity contact region 
formed between two spherical surfaces in adhesive con- 
tact. Greenwood and Johnson 30] have presented an el- 
egant approach, which I use in Appendix B to calculate 
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FIG. 17: The stress distribution in an asperity contact region. 
The stress is tensile close to the boundary line of the contact 
area. An attractive wall-wall interaction occur also outside 
the area of real contact and this interaction determines the 
detachment stress aa and the work of adhesion 7. 

Note that with P(-(Ta, C) 0, Eq. (39) takes the form 

F'iO^-a^iOA'iO (40) 

Let us consider the simplest case where is independent 
of the magnification In this case (40) gives 



dC 



[F(C)+(TaA(C)]-0 



or 



F(C) + (TaA(C) = ((To + (Ta)Ao 



(41) 



where we have used that F{1) = (TqAo and A{1) — Aq. 
Thus, if we assume that A{Q « Aq, a,s is typically the 
case for elastically hard solids and large magnification, 
then 



F{C) « ((To + (Ta)Ao 



(42) 



i.e., the normal force in the contact area is the sum of 
the applied load ctqAq and the adhesion load a^AQ. This 
result was derived in R ef. M using a different approach, 
and has often been used[22|- Note, however, that in most 
practical situations (Ta will depend strongly on the mag- 
nification and the relation (42) is then no longer valid. 
One exception may be systems with a suitable chosen 
graded elasticity, see Sec. 8.5. 

8.2 Elastic energy 

When complete contact occur between the solids on all 
length scales, Ud = U^^ with[lo| 



AoE 



4(1 



d^q qC{q) 



4(1 -z/2)' 

where qi = Ci is the largest surface roughness wave 
vector. In Ref. [lOj we accounted for the fact that only 
partial contact occur at the interface by using 

C^ci = 4(f!_^^2) 2^gg j^^' dC C'CiqoOPiO (43) 

where P{C) — ^(C)/^o is the ratio between the area of 
real contact (projected on the xy-plane), A{(^), when the 
system is studied at the magnification and the nomi- 
nal contact area ^o- Here we wih use a (shghtly) more 
accurate expression for jeS determined as foUows. 
First note that 



where Po{a,() = Po{(^) is the stress distribution assum- 
ing complete contact at all magnifications. 

The factor P(C) in Eq. (43) arises from the assumption 
that the elastic energy which would be stored (per unit 
area) in the region which we denote as "detached area" , 
if this area would not detach, is the same as the average 
elastic energy per unit area assuming complete contact. 
However, in reality the stress and hence the elastic energy 
which would be stored in the vicinity of the "detached 
area", if no detachment would occur, is higher than the 
average (that's why the area detach). Since the stored 
elastic energy is a quadratic function of the surface stress, 
we can take into account this "enhancement" effect by 
replacing P(C) in (43) with 



JZ^^^^da (a-ao)2p(a,C) 



or 



((^-^o)2>? 



(44) 



where in the denominator (...) is calculated assuming 
complete contact. Note that P{C) takes into account both 
the decrease in the contact area with increasing magnifi- 
cation, and the fact that the detached regions are the 
high tensile stress regions where —a > crdO- Thus, 
P{Q < P{C) and the amount of stored elastic energy at 
the interface will be smaller when using the ^(O-function 
instead of the P(C)-function. The P{C) function can be 
written in a convenient way, involving only P{C) rather 
than the distribution function P{cr, ^), as follows. Let us 
write 



{{a - ao)')c = F{0 - 2aoG(C) + a'oPiC) 
Consider first the integral 

/oo 
daa'PiaX) 



(45) 



Using (17) and that P(crc(C), C) = we get 
F'{0= doa-f 



/(C) 



da a 



, d^P 



Performing two partial integrations and using (35) gives 



F'iO = /(C) 



rIP 

^c(C)^(-^c(C),C) + 2P(C) 



(46) 



Integrating (17) over a gives 

P'(C) = -/(C)|^(-^c(C),C) 

Substituting this in (46) gives 

p'(C)-a2(c)P'(C) + 2/(C)P(C) 

In a similar way one can show that 



G'(C)-^(a)c = -acP'(C) 



(47) 



Thus 



-{{a ~ ao)')c = F'iO " 2aoG"(C) + alP'iC) 



= ('Tc(C)+ao)^P'(C) + 2/(C)P(C) 
Integrating this equation gives 

{i<^-^o?)c = dC [(a,(C') + aorP'iC) + 2/(C')^(C')] 

(48) 

Now, assume that complete contact occur at the interface 
at all magnifications. In this case P(C) = 1 and (48) takes 
the form 



(('^-<^o)'>C-/<2/(C') 
Thus, using (48) and (49) gives 

J^dC'aUOP'iC) , J^dC fiOPiC) 



2/NC7(C') 



/NC7(C') 



(49) 



(50) 



8.3 Numerical results 

In many applications surface roughness occur over 
many decades in length scale, e.g., from cm to nm (7 
decades in length scale). In this case it is not possible to 




h/ho 



FIG. 18; (a) The relative contact area at the highest mag- 
nification A{C^-i)/Ao, and (b) the (normalized) macroscopic 
efi'ective interfacial energy per unit area 7eff(l)/A7 as a func- 
tion of h/ho- We show the results using both the "old" form 
of the elastic energy and the improved form where is 
replaced by -P(C)- For an elastic sohd with a fiat surface 
in adhesive contact (at vanishing applied pressure ctq = 0) 
with a hard, randomly rough self affine fractal surface with 
the root-mean-square roughness ho = 0.5 /^m and the frac- 
tal dimension Df — 2.2. The long and short distance cut-off 
wavevectors go = 10^ m~^ and qi = 10^° m~^, respectively. 
The elastic solid has the Young's modulus E — 10 MPa and 
the Poisson ratio 0.5. The interfacial energy per unit area 
A7 = 0.05 J/m^ 

perform (numerically) integrals over the magnification di- 
rectly, but the integrals are easily performed if one change 
the integration variable ( — exp(/i) so that 

where fj,i — In^i . This change of variable is very impor- 
tant not only for the adhesional problems we consider 
here but also in other problems involving many length 
scales, e.g., rubber friction|3| or crack propagation fsills^ l 
in rubber. 

We now present an application relevant to the tire-rim 
sealing problemj^. A tire makes contact with the wheel 
at the rim, where the rubber is squeezed against the steel 
with a pressure of order 0.5 — 1 MPa. The contact be- 




FIG. 19: (a) The relative contact area A(()/Ao and (b) 
the (normalized) effective interfacial energy per unit area 
7cff(C)/A7 as a function of magnification. For the nominal 
pressure ao ~ 0.5 MPa. For the same system as in Fig. 1181 

tween the rubber and the steel must be such that neg- 
ligible air can leak through the non-contact rubber-steel 
rim area. In addition, the rubber-rim friction must be so 
high that negligible slip occur even during strong accel- 
eration or breaking where the torque acting on the tire 
is highest. Since the rubber- rim friction originate from 
the area of real contact, this require a large fractional 
contact area A{(i)/Aq. 

Consider an elastic solid with a flat surface in adhesive 
contact with a hard, randomly rough self affine fractal 
surface with the root-mean-square roughness (rms) ho = 
hrms = 0.5 /im and the fractal dimension Df — 2.2. The 
long and short distance cut-off wavevectors for the power 
spectrum go — 10^ m^^ and qi = 10^'^ m"^, respectively. 
This gives a surface area A* which is about 38% larger 
than the projected area Aq. 

The elastic solid has the (low-frequency) Young's mod- 
ulus E = 10 MPa and the Poisson ratio 0.5, as it typical 
for tire rubber. The interfacial binding energy per unit 
area A7 = 0.05 J/m^, as is typical for rubber in contact 
with a steel surface. 

Fig in^ a) shows the relative contact area A{(^i)/Ao and 
(b) the (normalized) effective macroscopic interfacial en- 
ergy per unit area 7off(l)/A7 as a function of h/ho, 



vanishing applied pressure ctq = 0. Here h is the rms 
roughness used in the calculation, i.e., the actual power 
spectrum has been scaled by the factor (h/ho)"^- We show 
results using both the "old" form of the elastic energy and 
the improved form where P{C) is replaced by P{C)- Note 
that the theory predict 7off(l) — for h — ho, i.e., the 
theory predict that the macroscopic adhesion vanish. 

In Fig. ^1 we show (a) the relative contact area 
A{(^)/Ao and (b) the (normalized) effective interfacial en- 
ergy per unit area 7eff(C)/^7 as a function of magnifica- 
tion for h = hf), when the surfaces are squeezed together 
with the nominal pressure (Jq — 0.5 MPa. Note that the 
adhesion increases the area of real contact (at the high- 
est magnification) A{(i) from « O.OQAq to ~ OAAAq. 
The adhesion becomes important already at relative low 
magnification C ~ 10 corresponding to the length scale 
A = 27r/(Cgo) ~ 10 //m. The magnification at the point 
where adhesion becomes important can be estimated by 
the condition that the nominal pressure ctq is of similar 
magnitude as the detachment stress: (Tq ~ {EA^qY^"^ 
which gives q ~ ao/EA-f w 5 x 10^ m~^ ot C, = q/qo ^ 5. 

Fig. [5D| shows the relative contact area (at the highest 
magnification) A{(i)/Ao as a function of the (nominal) 
squeezing pressure ctq with and without the adhesional 
interaction. For all squeezing pressures, the area of real 
contact is much larger when the adhesional interaction is 
included. 



0.4 



A/Af 



0.2 



with adhesion 




0.1 0.2 0.3 
pressure (MPa) 



0.4 



0.5 



FIG. 20: The relative contact area (at the highest magni- 
fication) A((i)/Ao as a function of the (nominal) squeezing 
pressure ctq with and without the adhesional interaction. For 
the same system as in Fig. 1181 



Surface roughness is an important factor which influ- 
ence the rate of leakage through seals. If one assume 
that the contact regions at any magnification are ap- 
proximately randomly distributed in the apparent con- 
tact area, we expect from percolation theory that the 
non-contact region will percolate when A/Aq ~ 1 — Pc, 
where pc ~ 0.7 is the site percolation number. Since 
(because of adhesion) at the highest (nanometer) resolu- 



tion, A/Ao « 0.44 > 1 — pc ~ 0.3, we conclude that no 
roughness induced gas leakage will occur in the present 
case. 

Next, let us focus on the tire-rim friction. The nom- 
inal contact area between the tire rim and the steel is 
roughly 0.05 m^. The kinetic frictional stress when rub- 
ber is sliding on a perfectly smooth clean substrate de- 
pend on the slip velocity but is typically of order mm 
« 0.25 MPa. Thus, if the rubber is in perfect contact 
with the steel, the shear force « 1.3 x lO'' N. For the 
system studied above we found A/Aq « 0.44 giving the 
friction force « 5.7 kN. If the load on a tire is 3.5 kN 
and the tire-road friction coefficient of order unity, the 
maximal friction force -Ff ~ 3.5 kN which is below the 
estimated rubber-steel friction force, indicating negligi- 
ble slip of the tire relative to the steel. However, slip 
is sometimes observed, indicating that the assumptions 
made above about, e.g., the steel surface roughness, may 
differ from real systems. It is clear, however, that analysis 
of the type presented above, using the measured surface 
topography for the steel (or aluminum) surface, and the 
measured rubber viscoelastic modulus E and interfacial 
rubber-steel binding energy per unit area A7, may help 
in the design of systems where no (or negligible) tire slip 
(and gas leak) occur at the rim. 

8.4 Detachment stress 

The expression (36) for the detachment stress will, in 
general, break down at very short length scales (high 
magnification). This can be understood as follows. For 
"simple" solids the elastic modulus E can be related to 
the surface energy (per unit area) 71 ~ Ea, where a is 
an atomic distance (or lattice constant). Thus (36) pre- 
dict the that detachment stress at the highest (atomic) 
resolution a, where jcS = A7 and 5 ^ 1 /a, is of order 

(A77i)i/Va 

However, at the atomic scale the detachment stress must 
be of order a* = A'j/a. For many elastically stiff solids 
with passivated surfaces 71 >> A7. For such cases we 
suggest to use as an interpolation formula 



where 



f^a(C) 



f a7off(C)^^A 



1/2 



We may define a cross-over length 1^ via the condition |3£ 

1/2 



c?E 
A7 



For elastically stiff solids with passivated surfaces, the 
the cross-over length Ic may be rather large, e.g., of or- 
der ~ 1 for silicon with the surfaces passivated 
by grafted layers of chain molecules. However, for elas- 
tically soft solids such as the PMMA adhering to most 
surfaces (see Sec. 10.1) A7 ^ ^ 0.1 J/m^ so that 
the correction above becomes negligible. In this latter 
case, with a ~ 1 A the detachment stress at the highest 
(atomic) magnification cra(Ci) ~ 1 GPa. More gener- 
ally, for unpassivated surfaces (even for stiff materials) 
71 « A7 and the correction becomes unimportant. How- 
ever, most surfaces of hard solids have passivated surfaces 
under normal atmospheric conditions and in these cases 
the "correction" above may become very important. 

8.5 Flaw tolerant adhesion 

The breaking of the adhesive contact between solids 
is usually due to the growth of interfacial cracks which 
almost always exist due to the surface roughness (see 
above) or due to other imperfections, e.g., trapped dust 
particles. Since the stress at a crack tip may be much 
larger than the nominal (or average) stress at the inter- 
face during pull-off, the adhesive bond will usually break 
at much lower nominal pull-off stress than the ideal stress 
a* — o'a(Ci) ~ ^7/^1 where a is a typical bond distance 
of order Angstrom. An exception is the adhesion of very 
small solids (linear size I < Ic), where the pull-off force 
ma y a p proach the ideal value a*, see Sec. 8.4 and Ref. 

It has been suggested that the strength of adhesion 
for macroscopic solids can approach the ideal value a* 
if the material has an elastic modulus which increases 
as a function of the distance z from the surface (graded 
elasticity') (s?! . This can be understood as follows: The 
detachment stress for an interfacial crack of linear size A 



A(l-Z/2) 



1/2 



(51) 



Since the strain field of a crack of linear size A extend 
into the solid a distance ~ A, it will experience the elas- 
tic modulus for z w A. Thus, the ratio E/\ in (51) will 
be independent of the magnification ~ L/X if E{z) ^ z. 
Hence if ^^s would be independent of the magnification 
C, then the detachment stress (Ta(C) would also be inde- 
pendent of the magnification and of order the theoretical 
maximum a*. This argument was presented by Yao and 
Gaols^- However, our calculations show that in most 
cases 7cff will increase rapidly with increasing magnifi- 
cation, in which case the elastic modulus must increase 
faster than z with the distance from the surface in order 
for Ua to stay constant. In this case there will be no "uni- 
versal" optimal form of the z-dependence of E{z), since 
it will depend on the surface roughness power spectrum. 

8.6 Contact mechanics and adhesion for vis- 
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FIG. 21: Asperity contact region during pull-off of a rubber 
block from a rough surface. The dashed area indicate the 
dissipative region at the opening crack. The theory of crack 
propagation in viscoelastic solids is only valid as long as the 
linear size d of the dissipative region at the crack tip is smaller 
than the linear size D of the asperity contact area. 



coelastic solids 

The contact between viscoelastic solids and hard, ran- 
dom ly r ough, substrates has been studied in Ref. 
and jl^ for stationary contact, and in Ref. for slid- 
ing contact. In these papers the adhesional interaction 
between the solids was neglected. No complete theory 
has been developed for viscoelastic contact mechanics be- 
tween randomly rough surfaces including the adhesion. 
However, some general statements can be made. 

Assume that the loading or pull-off occur so slowly 
that the viscoelastic solid can be considered as purely 
elastic everywhere except close to the asperity contact 
regions. That is, we assume that in the main part of the 
block the deformation rate is so slow that the viscoelastic 
modulus E{uj) « £^(0). However, the detachment lines 
can be considered as interfacial cracks and close to the 
crack edges the deformation rate may be so high and the 
full viscoelastic modulus must be used in these regions of 
space. 

The adhesional interaction is in general much more 
important during separation of the solids than during 
the formation of contact (loading cycle). This is related 
to the fact that separation involves the propagation of 
interfacial opening cracks, while the formation of con- 
tact (loading cycle) involves the propagation of interfacial 
closing cracks. For flat surfaces, the effective interfacial 
energy (per unit area) necessary to propagate opening 
cracks is of the form 13 iL 

-f{v)^Aj[l + fiv)] 

where f{v) is due to the viscoelastic energy dissipation in 
the vicinity of the crack tip, and depend on the crack tip 
velocity v. For rubber-like viscoelastic solids, f(v) may 
increase with a factor of ^ 10"^ with increasing crack- 
tip velocity. Thus, the viscoelastic deformations of the 
solid in the vicinity of the crack tip may strongly in- 



crease the pull-off force, as compared to the purely elas- 
tic case where 7 = A7. We note, however, that for a 
system where the contact only occur in small (randomly 
distributed) asperity contact regions, the formula above 
will only hold as long as the linear size of the dissipative 
region is smaller than the diameter of the contact area, 
see Fig. EH 

The effective interfacial energy gained during the prot 
agation of closing cracks is of the approximate formji 



7(1;) 



A7 



1 + 



Thus, in this case "f{v) is always smaller than (or equal 
to) A7, so that the adhesional interaction is much less 
important than during pull-off. 
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FIG. 22: Terminal elements (circles) in animals with hairy 
design of attachment pads. Note that heavier animals exhibit 
finer adhesion structures. Adapted from Ref. |4ll| . 



8.7 Strong adhesion for (elastically) stiff mate- 
rials: Natures solution 

In order for an clastic solid to adhere to a rigid rough 
substrate, the elastic modulus of the solid must be small 
enough. The physical reason is that the surface of the 
elastic solid must deform elastically for contact to occur 
at the interface, and if the elastic energy necessary for the 
deformation is larger than the interfacial binding energy 
A7A1 (where Ai is the area of real contact) no adhesion 
will be observed. Thus, experiment and theoretical stud- 
ies have shown that already a surface roughness with a 
root-mean-square amplitude of order a few micrometer 
is enough in order to kill adhesion even for very soft ma- 
terials such as rubber with a (low frequency) modulus of 
order ~ 1 MPa. In spite of this strong adhesion is possi- 
ble even for intrinsically relative stiff materials if instead 
of compact solids, the solids (or thick enough layers at the 
surfaces of the solids) form foam-like or fiber-like struc- 
tures; the effective elastic properties of such non-compact 
solids can be much smaller than for the compact solids. 



This fact is made use of in many biological adhesive sys- 
tems. 

The adhesive microstructures of lizards and many in- 
sects is the result of perhaps millions of years of devel- 
opment driven by the principle of natural selection. This 
has resulted in two different strategies for strong adhe- 
sion: foam-like structures as for the adhesion pads of 
cicada (see Fig. I14|l . and fibers or hierarchical small- 
fibers-on-large-fibers structures (usually terminated by 
thin elastic plate-like structures), as for lizards, flies, bugs 
and spiders, see Fig. 123 These adhesive systems are 
made from keratin-like proteins with an elastic modulus 
of order ~ 1 GPa, i.e., ~ 1000 times higher than for rub- 
ber, but the effective elastic modulus of the non-compact 
solids can be very low^^, making strong adhesion pos- 
sible even to very rough surfaces such as natural stone 
surfaces. 

9 Elastoplastic contact mechanics with adhesion 

It is possible to apply the contact mechanics theory to 
the rather complex but very important 42J case of elasto- 
plastic contact with adhesion. In this case one need to 
solve the equation for the stress distribution 

with the boundary conditions derived above 



^'(^Y,C) = 0, 

and the "initial" condition 



P(-(Ta,C)=0 



P{<J, 1) = S{<J - do) 

If one introduce the new variable 

<7y — a 
s = 

and consider P(cr, C) = P{s,Q as a function of s and (, 
and denote crb(C) = +o'a(C)i then the equations above 
takes the form 



dP _a'^dP ^ fipd^P 
9C o'b ds ds^ 



p(o,c) = 0, 



p(i,C) = o 



(52a) 



(526) 



P(s,l) =J(aY-fTo-sab(l)) (52c) 
The boundary conditions (52b) are satisfied if one expand 



P(s,C) = ^A„(C) sin( 



Trns 



Substituting this in (52a) gives a linear system of first 
order ordinary differential equations for An (C) which can 
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FIG. 23: The calculated relative contact area as a function 
of the (nominal) pressure ctq including adhesion (upper solid 
hne) and without adhesion. Squares: the normalized, (nom- 
inal) shear stress r/rf for PMMA sliding on TMS as a func- 
tion of the (nominal) pressure. The nominal shear stress r 
has been divided by T{ = 50 MPa, which is the (measured) 
shear stress in the area of real contact (see text). Circles: the 
normalized, (nominal) shear stress r/rf (with T{ — 5 MPa) 
for PMMA sliding on OTS as a function of the (nominal) 
pressure. 

be integrated and solved using matrix inversion. We will 
present the detailed derivation and numerical results else- 
where. 

10 Applications 

I present two applications involving very smooth sur- 
faces, prepared by cooling liquid glassy materials be- 
low the glass transition temperature (see Sec. 2). We 
first consider the contact between a PMMA lens and a 
flat hard substrate, and discuss the experimental data 
of Bureau et al0|- also consider contact mechan- 
ics for Pyrex glass surfaces, where plastic yield occur in 
the asperity contact regions, and discuss the results of 
Restagno et al|4^. 

10.1 Contact mechanics for PMMA 

In the experiment by Bureau et al.(43|], a PMMA lens 
was prepared by cooling a liquid drop of PMMA from 
250 °C to room temperature. In the liquid state the sur- 
face fluctuations of vertical displacement are caused by 
thermally excited capillary waves (ripplons). When the 
liquid is near the glass transition temperature Tg (about 
100 °C for PMMA), these fluctuations become very slow 
and are finally frozen in at the glass transition. Thus, 
the temperature T in (1) is not the temperature where 
the experiment was performed (room temperature), but 
rather the glass transition temperature Tg. 

The (upper) solid line in Fig. 1231 shows the calculated 
relative contact area as a function of the (nominal) pres- 
sure ctq = F^/Aq (where is the squeezing force). In 



the calculation we have used the equations above with 
the measured elastic (Young) modulus E = 2.9 GPa 
and surface tension 7 = 0.04 J/m^. We have also used 
the glass transition temperature Tg « 370 K, the short- 
distance cut-off wavevector qi ~ 7 x 10^ m"-'^ and Qc = 
1.7 X 10^ m^^ (see Eq. (1)). The cross-over wavevector 
Qc (or the bending stiffness k) has not been measured for 
PMMA, but has been measured for other systems. Thus, 
for alkanes at T « 100 °C for C20 and C36, ~ 4.4 x 10^ 
and « 2.7 x 10^ m"-'^, respectively [isf . There are also 
some studies of polymers [s^ 0| (using the atomic force 
microscope) giving gc ~ 1 x 10® m^^. The interfacial 
binding energy A7 can be estimated using 20] A7 « 
2(7172)^/^ ~ 0.07 J/m^, where 71 = 7 « 0.04 J/m^ is the 
surface energy of PMMA and 72 ~ 0.02 J/m^ the surface 
energy of the (passivated) substrate. In the calculations 
we used A7 = 0.06 J/m^. 

If one assumes that, in the relevant pressure range (see 
below) , the frictional shear stress T{ between PMMA and 
the substrate is essentially independent of the normal 
stress a in the asperity contact regions, then the calcu- 
lated curve in Fig. |531is also the ratio between the nom- 
inal (or apparent) shear stress r and the true stress rf; 
A/Aq = t/t{ (note: the friction force F{ — TfA — tAq). 

The circles and squares in Fig. [221 show the measured 
data of Bureau et al.ji^l for r/rf. They performed ex- 
periments where a PMMA lens, prepared by cooling a 
liquid drop of PMMA, was slid on silicon wafers (which 
are nearly atomically smooth) covered by a grafted silane 
layer. Two different types of alkylsilanes were employed 
for surface modification, namely a trimethylsilane (TMS) 
and octadecyltrichlorosilane (OTS). The squares in Fig. 
123 are the shear stress, divided by rf = 50 MPa, for 
PMMA shding on TMS as a function of the (nominal) 
pressure ctq- The shear stress rf = 50 MPa was deduced 
from multi-contact experiments [43 (where the local pres- 
sure in the asperity regions is so high as to give rise 
to local plastic deformation) using rf « rH/ao, where 
H sa 300 MPa is the hardness (yield stress) of PMMA as 
determined from indentation experiments. This equation 
follows from tAq = TfA and AH = Aoao. The circles in 
Fig. 031 show r/rf for PMMA sliding on OTS as a func- 
tion of the (nominal) pressure ctq. For this system no 
multicontact experiments were performed, and we have 
divided the apparent shear stress by rf = 5 MPa, chosen 
so as to obtain the best agreement between the theory 
and experiment. Hence, in this case rf represent a theo- 
retically predicted shear stress for the multicontact case; 
we suggest that multicontact experiments for the system 
PMMA/OTS are performed in order to test this predic- 
tion. 

The calculation presented above shows that the lo- 
cal pressure in the asperity contact regions is of or- 
der ~ 2 GPa, which is much higher than the macro- 
scopic yield stress (about 300 MPa) of PMMA. However, 
nanoscale indentation experiments [48i| have shown that 
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FIG. 24: The adhesion force (or pull-off force) Fad as a func- 
tion of the squeeze force _Fn to the power 1/3. The strait line 
is a fit to the experimental data. Adapted from Ref. |44j 

on the length scale (~ 10 nm) and indentation depth 
scale (~ 1 nm) which interest us here, the indentation 
hardness of PMMA is of order ~ 10 GPa, and no plastic 
yielding is expected to occur in our application. 

The study above has assumed that the frictional shear 
stress at the PMMA-substrate interface is independent 
of the local pressure p = cr, which varies in the range 
p « 50 — 100 MPa. This result is expected because the 
shear stress will in general exhibit a negligible pressure 
dependence as long as p is much smaller than the ad- 
hesional pressure Pad- The adhesional pressure is de- 
fined as follows: In order for local slip to occur at the 
interface, the interfacial molecules must pass over an 
energetic barrier of magnitude 6e. At the same time 
the spacing between the surfaces must locally increase 
with a small amount (some fraction of an atomic dis- 
tance) Sh. This correspond to a pressure work paP'Sh, 
where is the contact area between the molecule, or 
molecular segment, and the substrate. Thus the total 
barrier is of order (5e + pa^6h = a^Sh^pi^^d + p) where 
Pad — Se/{a^Sh). For weakly interacting systems one 
typically have [ssl Se/a'^ w 1 meV/A and Sh « 0.01 A 
giving Pad ~ 10 GPa, which is at least one order of magni- 
tude larger than the pressures in the present experiment. 
We note that the pressure Pad is similar to the pressure p 
which must be applied before the viscosity of liquid hy- 
drocarbon oils start to depend on the applied pressure. 
The reason for the pressure dependence of the viscosity 
of bulk liquids is similar to the problem above, involving 
the formation of some local "free volume" (and the as- 
sociated work against the applied pressure) in order for 
the molecules to be able to rearrange during shear. 

10.2 Contact mechanics for Pyrex glass 

Cohesion effects are of prime importance in powders 
and granular media, and they are strongly affected by 



the roughness of the grain surfaces. For hydrophilic and 
elastically stiff solids such as silica glass, the cohesion 
is usually strongest under humid conditions because of 
the capillary forces exerted by small liquid bridges which 
form at the contact between grains. Recently Restagno et 
al0 have shown that adhesive forces can be observed 
for stiff solids even in the absence of capillary bridges 
if the surfaces are smooth enough. However, even for 
very smooth Pyrex glass surfaces (see below) the adhe- 
sional force between a small glass ball and a flat glass 
surface is a factor ~ 10"'^ smaller than expected for per- 
fectly smooth surfaces. The elastic modulus of Pyrex 
glass is about ~ 10^ times higher than for PMMA, and 
for this case the theory (see below) predict that even an 
extremely small surface roughness will result in no adhe- 
sion, assuming that only elastic deformation occur in the 
contact areas between the Pyrex glass surfaces. However, 
the high contact pressure in the asperity contact regions 
will result in local plastic yielding, and if the surfaces are 
squeezed together with a high enough force, the plasti- 
cally flattened regions will be so large that the adhesional 
interaction from these contact regions will overcome the 
(repulsive) force arising from elastic deformation of the 
solids in the vicinity of the contact regions, resulting in 
a non- vanishing pull-off force. This effect was observed 
in the study by Restagno et al.0|- 

In the experiments by Restagno et al.0| a sphere with 
radius R = 1.8 mm was squeezed against a flat substrate, 
both of fire-polished Pyrex. The surface topography was 
studied using the AFM over an area 10 x 10 fim^. The 
maximum height fluctuation was of order ~ 1 nm as 
expected from the frozen capillary waves. The adhe- 
sion experiments where performed in a liquid medium 
(dodecane) in order to avoid the formation of capillary 
bridges. If the ball and flat was squeezed together by 
a force F < Fc — 10^® N and then separated, no ad- 
hesion could be detected. However, if the squeeze force 
was larger than this (critical) value, a flnite pull-off force 
could be detected. 

Fig. shows the measured pull-off or adhesion force 

1/3 

Fad as a function of F-^ , where Fn is the force with 
which the ball and the flat was squeezed together (i.e., 
the maximum force during the load cycle). Note that 

1 /3 

Fad depends approximately linearly on F-^ . As shown 
by Restagno et al.0| this can be simply explained us- 
ing standard contact mechanics for a ball-on-flat conflg- 
uration, and assuming that plastic yielding occur in all 
asperity contact regions. 

For elastically stiff solids such as Pyrex, the pull-off 
force for a perfectly smooth spherical object in adhesive 
contact with a perfectly smooth flat substrate is given 
approximately by the Derjaguin approximation 21] : 

Fad - 2ttAjR (53) 

where A7 = 27s; is twice the interfacial energy (per unit 
area) of the Pyrex-dodecane system. The same equation 



is valid when surface roughness occur, assuming that the 
dominant roughness has wavelength components smaller 
than the nominal contact area between the solids. How- 
ever, in this case 7 must be replaced with "/eS, which is 
an effective interfacial energy which takes into account 
the surface roughness. The roughness modifies 7 in two 
ways. First, the area of real contact A is smaller than 
the nominal contact area Aq. This effect was taken into 
account by Restagno et al. by writing the (macroscopic, 
i.e., C = 1) Toff = 7^Mo- The second effect is the elastic 
energy stored in the vicinity of the contacting asperi- 
ties. This elastic energy is given back during pull-off and 
will therefore reduce 7off. This effect was neglected by 
Restagno et al. 

During squeezing the nominal contact area is given 
(approximately) by Hertz's theory 
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where E* — 2E/3{1 — v'^). Restagno et al. assumed that 
plastic yielding has occurred in the area of real contact 
so that A — F/H where H is the indentation hardness 
of the solids. Thus using (54): 
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Combining this with (53) and (55) gives 
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The strait line in Fig. 1241 is given by (56) with 7 = 7si = 
0.43 J/m^, and with R — 1.8 mm, and the known elastic 
modulus, Poisson ratio and Pyrex hardness. 

The theory above is consistent with the experimen- 
tal data but does not address the following fundamental 
questions: 

(a) Why is no adhesion observed when the solids are 
brought together "gently" , i.e., in the absence of the load- 
ing cycle (squeeze force F^ ~ 0)? 

(b) Why does adhesion occur when the solids are first 
squeezed into contact. That is, why does not the stored 
elastic energy at the interface break the interfacial bonds 
during removal of the squeezing force F^l 

We will now address these questions using the contact 
mechanics theory for randomly rough surfaces described 
above. Assume that the surface roughness on the Pyrex 
surfaces is due to frozen capillary waves. In fig. I^S] (up- 
per curve) we show the surface roughness power spectrum 
calculated using (1) with Tg = 1100 K, 7 = 1.0 J/m^ and 
= 0.5 X 10^° m~^. The lower curves in Fig. I^BT a') 
and (b) show the relative contact area A/Aq (at the 
highest magnification C = Ci) ^.nd the effective macro- 
scopic (i.e., for C = 1) interfacial binding energy (per 
unit area) 7off(l) as a function of the root-mean-square 
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FIG. 25: The surface roughness power spectrum of the orig- 
inal (upper curve) and the plastically smoothened surface. 
For Pyrex glass using the power spectrum calculated from 
the model with frozen capillary waves. In the calculation of 
the smoothened power spectrum we used the elastic modu- 
lus E = 32 GPa and Poisson ratio v = 0.2 and the hardness 
CTY = 6 GPa. 



surface roughness h divided by Hq. Here Hq is the root- 
mean-square roughness of the actual surface and h the 
rms roughness used in the calculation. Thus, the power 
spectrum has been multiplied (or scaled) with the factor 
(h/ho)'^ . The macroscopic interfacial energy 7off(l) in (b) 
is in units of the bare interfacial energy A7. Note that 
7cft(l) and the contact area both vanish for h > 0.94/io, 
so that no adhesion is expected in this case. This is in 
accordance with the experiment since no adhesion is ob- 
served if the glass surfaces are brought together without 
squeezing. 

Fig. 1^ shows the calculated relative contact area as 
a function of the magnification ^ using C(g) from Fig. 

as the surfaces are squeezed together with the nom- 
inal pressure ctq — 20 MPa. In (a) we show the clas- 
tic contact area both for the elasto-plastic model with 
the hardness ay — & GPa (lower curve), and for infinite 
hardness, where only elastic deformation occur (upper 
curve). In (b) we show the plastic contact area. Note 
that the magnification ^ — 10'' correspond to the wave 
vector q — (qo — 2n/X — 10^° m~^, i.e., to atomic 
resolution A « 6 A. The plastic yielding occur in a 
relative narrow range of resolution; for C w 3000, cor- 
responding to A « 2 nm, about half the contact area 
has yielded. Thus a plastic yielded solid-solid asper- 
ity contact region has a typical area of order 4 nm^. 
Note that in the experiment (Fig. I24|l the nominal pres- 
sure .Fn /Aq 



F^^^{E* /R)'^/^ /tt, so the nominal pressure 
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20 MPa correspond to F^^^ = 0.06 N^/^. 



Fig. 1271 shows that the short-ranged surface roughness 
will be plastically fiattened so that after squeezing the 
surfaces will conform at the shortest length-scale. We 



A/A( 



0.5 





^ (a). 




plastically^ 




smoothened 




^surface 


original 




surface 





Teff /Ay 




FIG. 26: The calculated relative contact area A/Ao (at the 
highest magnification) and the effective interfacial binding en- 
ergy (per unit area) j^s as a function of the surface rough- 
ness h/ho- Here ho is the root- mean-square roughness of the 
actual surface and h the rms roughness used in the calcula- 
tion. Thus, the power spectrum has been multiplied with the 
factor {h/hoY . (a) The contact area area both for the plas- 
tically smoothened surface, and for the original surface, (b) 
The macroscopic interfacial energy 7cff(l) in units of the bare 
interfacial energy A7. For Pyrex glass using the power spec- 
trum calculated from the model with frozen capillary waves. 
In the calculation we used the elastic modulus E = 32 GPa 
and Poisson ratio u = 0.2 and A7 = 0.43 J/m^. 



can take this into account in the analysis of the pull-off 
force by replacing the power spectrum C (q) with a power 
spectrum C{q) corresponding to a "smoothened" surface, 
where the short wavelength roughness is removed. One 
approach would be to cut-off C (g) for q > q^, where ql — 
goCi where Q is the magnification where, say 50% of the 
contact area has yielded plastically. Another approach is 
to define 

C{qoC) =(^l- ^) C{qoO 

where A^^ = Fn / cty is the contact area when the whole 
contact area has yielded plastically. In Fig. [53 (lower 
curve) we show C{q) using this definition. In Fig. [^(a) 
and (b) (upper curves) we show the the relative contact 
area A{(i)/Ao and the macroscopic effective interfacial 



o 
< 

CO -4 
o 



(a) 






no plastic" 




yield , 








FIG. 27: The calculated relative contact area as a function 
of the magnification (. (a) The elastic contact area both for 
the elasto-plastic model with the hardness cty = 6 GPa (lower 
curve), and for infinite hardness, where only elastic deforma- 
tion occur (upper curve), (b) The plastic contact area. For 
Pyrex glass using the power spectrum calculated from the 
model with frozen capillary waves. In the calculation we used 
the elastic modulus E/2 = 32 GPa (the factor 1/2 arises from 
the fact that both solids are Pyrex with the same elastic mod- 
ulus E) and Poisson ratio u = 0.2. The magnification ^ = 10* 
correspond to the wave vector q = (^qo = 27r/A = 10^° m"^, 
i.e., to atomic resolution A ~ 6 A. The plastic yielding oc- 
cur in a relative narrow range of resolution; for ( ~ 3000, 
corresponding to A « 2 nm, about half the contact area has 
yielded. Thus a plastic yielded solid-solid contact region have 
typical area of order 4 nm^. In the calculation the nominal 
squeezing pressure ctq = 20 MPa. 



energy (per unit area) 7ofr(l) as a function of h/ho for 
the smoothened surface. Note that in this case for h = ho 
there is nearly perfect contact at the interface, and the 
macroscopic 7off is about half of A7. In the actual ap- 
plications the short- wavelength roughness is not removed 
(or smoothened) over the whole surface, but only in small 
macro-asperity contact regions, with a total area propor- 
tional to the maximum squeezing force i^N during the 
load cycle. Nevertheless, the calculation shows that one 
may expect nearly perfect contact within the smoothened 
area, and only a relative small reduction in 7eff due to 
the stored elastic energy at the interface. In fact, the 
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FIG. 28: An elastic solid (dotted area) in adhesive contact 
with a hard substrate (black) with roughness on two diflterent 
length scales, (a) The (short-wavelength) substrate roughness 
is so high that no contact occur and the effective interfacial 
energy y^s = 0. (b) The small-wavelength roughness is flat- 
tened by plastic yield close to the top of the large asperity by 
an earlier squeezing cycle (loading force Jn). The contact area 
at zero load is determined by the adhesive interaction which 
is able to pull the surfaces into complete contact within the 
the smoothened region with an area ~ _Fn. The contact area 
cannot grow beyond the smoothened region because of the 
short wavelength roughness outside the smoothened region- 
we refer to this as the roughness barrier effect, (c) In the 
model calculation the surface is smoothened everywhere and 
the solids make contact nearly everywhere. 

reduction in 7off(l) will be smaller than indicated in Fig. 
|2ni because the elastic deformations will be smaller when 
the solids only makes contact close to the (smoothened) 
macro-asperity tops (compare Fig. |5Sl(b) and (c)). 

In Fig. [2H1 we illustrate (schematically) the basic re- 
sults of our study. Here an elastic solid (dotted area) 
in adhesive contact with a hard substrate (black) with 
roughness on two different length scales. For the origi- 
nal surface (not plastically smoothened) (a) the (short- 
wavelength) substrate roughness is so high that no con- 
tact occur and the effective interfacial energy 7cff — 0. 
During the loading cycle, the short wavelength rough- 
ness at the top of the macro-asperities is flattened by 
plastic yield. For the smoothened surface the contact 
area at zero load is determined by the adhesive inter- 
action which is able to pull the surfaces into complete 
contact within the the smoothened region (see (b)), with 



an area proportional to the maximum loading force Fn 
during the loading cycle. The contact area cannot grow 
beyond the smoothened region because of the short wave- 
length roughness outside the smoothened region-we refer 
to this as the roughness barrier effect. In the model cal- 
culation the surface is smoothened everywhere as in (c) 
and the solids makes contact nearly everywhere. 

11 On the philosophy of contact mechanics 

Recently it has been argued that contact mechanics 
models based on continuum mechanics cannot be used 
to determine the interfacial pressure distribution with 
atomistic resolution, or the adhesion or friction between 
the contacting solids, since these properties depend on 
atomic scale surface roughness (surface adatoms, steps, 
and so on) . But the aim of continuum mechanics was 
never to describe the physics at the atomistic level but 
only at larger length scales. Thus, a fundamental ap- 
proach to the problems mentioned above consist of ap- 
plying continuum mechanics down to a length scale of 
order a few nanometer. At this length scale or resolu- 
tion continuum mechanics will correctly predict the con- 
tact area and the pressure distribution. At shorter length 
scale other methods such as molecular dynamics or quan- 
tum mechanical methods must be employed which takes 
into account the atomistic and chemical nature of real 
surfaces. However, in most practical applications such 
a detailed approach will not be possible because too lit- 
tle is usually known about the actual interface, e.g., thin 
(nanometer or less) contamination layers may dominate 
the physics at the atomistic level. Nevertheless, even in 
these cases the continuum mechanics approach may be 
extremely useful as it determines the size and the pres- 
sure distribution in the effective contact areas where most 
of the interesting physics takes place, see Fig. I^Hl 

12 Comment on numerical studies of contact 
mechanics 

Many numerical studies based on the finite element 
method, or other numerical methods, have been pre- 
sented for the elastic or elastoplastic contact between 
randomly rough (e.g., self affine fractal) surfaces |25l Isol 
l^ls^. However, numerical methods cannot be used to 
study the influence of roughness on the contact mechan- 
ics for macroscopic systems involving perhaps 7 decades 
in length scales (from nanometer to cm). The compu- 
tational time, and the required memory space, scales at 
least linearly with the number of grid points N — {L/af' 
in the a;?/-plane (a is the grid size and L the linear size 
of the system studied). (The linear scaling ^ N assumes 
using a multi-scale approach'53].) Thus, increasing the 
linear size of the system with one decade L lOL will, 
at least, increase the needed memory size and the compu- 
tational time by a factor of 100, and at present converged 
(see below) numerical studies are limited to systems with 




FIG. 29: A fundamental approach to contact mechanics, ad- 
hesion and friction consist of applying continuum mechanics 
down to a length scale of order a few nanometer. At this 
length scale or resolution continuum mechanics will correctly 
predict the contact area and the pressure distribution. At 
shorter length scale other methods such as molecular dynam- 
ics or quantum mechanical methods must be employed which 
takes into account the atomistic and chemical nature of real 
surfaces. 
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FIG. 30; The exact and approximate stress distribution 
P{a,Q. The exact P{a,Q vanish for a = 0, but this re- 
sult is not observed in numerical calculations owing to a too 
sparse integration grid (see text for details). The two stress 
distribution functions correspond to the same load, but the 
approximate stress distribution gives 22% larger area of con- 
tact. 

2 or maximum 3 decades in length scale. Nevertheless, 
numerical studies even on relative small systems may be 
useful to test the accuracy of analytical contact mechan- 
ics theories. 

All numerical studies of contact mechanics for ran- 
domly rough surfaces I am aware of are non-converged 
in the following sense. Assume that the wavelength of 



the shortest surface roughness component is Ai. In order 
to properly describe contact mechanics it is necessary to 
have at least 10 x 10 grid points and preferably much 
more (maybe 100 x 100) within one asperity contact re- 
gion which may be smaller than, or of order Ai x Ai. 
Only when this smallest Hertzian-like contact problem is 
accurately described can one obtain accurate results for 
the area of real contact. It is easy to show that when too 
few grid points are used, the numerical study will give an 
overestimation of the area of real contact. The reason for 
this is that the true stress distribution vanishes at zero 
stress (see Sec. 4), and in order to obtain this result (ap- 
proximately) in numerical studies it is necessary to have 
such a high density of grid points that the contact region 
close to the non-contact area is accurately described (it 
is this region which determine the boundary condition 
P(0, C) = 0). In all reported studies I am aware of this is 
not the case and the stress distribution remain finite (and 
large) at cr = as indicated in Fig. 1201 However, since 
the applied load is determined by the first moment of 
the pressure distribution, 

/•oo 
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it follows that the numerically determined P{(7, C) curve 
must be below the exact stress distribution curve for 
larger a in order for i^N to be the same. But since the 
area of real contact A is given by 
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it follows that the numerically determined contact area 
will be larger than the actual one. 

13 Summary and conclusion 

I have reviewed basic contact mechanics theories. For 
surfaces with random roughness on many different length 
scales, the contact area (in the absence of adhesion) is 
proportional to the squeezing force at low enough squeez- 
ing force. In addition, theory and numerical simulations 
show that when the contact area is proportional to the 
squeezing force, the stress distribution at the interface 
is independent of the squeezing force and the distribu- 
tion of sizes of the contact regions does not depend on 
the squeezing force. Thus, when the squeezing force in- 
creases, new contact regions are formed in such a way 
that the distribution of contact regions and the pressure 
distribution remains unchanged. This is the physical ori- 
gin of Coulomb's friction law which states that the fric- 
tion force is proportional to the normal (or squeezing) 
force 2] , and which usually holds accurately as long as the 
block-substrate adhesional interaction can be neglected. 

I have presented a detail discussion about the bound- 
ary conditions obeyed by the stress probability distribu- 
tion P{a, C) for elastic, elastoplastic and adhesive contact 



between solids. These boundary conditions are needed 
when studying contact mechanics using the formaUsm 
developed in Ref. 0- I have presented numerical re- 
sults illustrating some aspects of the theory. I have ana- 
lyzed contact mechanics problems for very smooth poly- 
mer (PMMA) and Pyrex glass surfaces prepared by cool- 
ing liquids of glassy materials from above the glass tran- 
sition temperature, and shown that the surface rough- 
ness which results from capillary waves when a glassy 
material is cooled below the glass transition tempera- 
ture can have a large influence on the contact mechan- 
ics between the solids. The analysis suggest a new ex- 
planation for puzzling experimental results [L. Bureau, 
T. Baumberger and C. Caroli, arXiv:cond-mat/0510232 
about the dependence of the frictional shear stress on 
the load for contact between a glassy polymer lens and 
flat substrates. I have also studied the contact mechanics 
between Pyrex glass surfaces where plastic deformation 
occur at the nanometer length scale. Finally, I have given 
some comments about the possibility of testing analyti- 
cal contact mechanics theories using numerical methods, 
e.g., fine element calculations. 

Acknowledgement I thank T. Baumberger and L. 
Bureau for drawing my attention to their interesting ex- 
perimental data and for supplying me the data used in 
Fig. I23| and E. Charlaix for drawing my attention to 
Ref. (4^]. I thank Pirelli Pncumatici for support. 

Appendix A: Power spectrum of frozen capil- 
lary waves 

Surfaces of liquids have thermally excited capillary 
waves. There are three contributions to the potential 
energy of capillary waves: a gravitational contribution, 
a surface tension contribution and another contribution 
which depend on the curvature of the liquid surface. The 
surface tension contribution depend on the surface area 
and can be written as 



If we write 
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The gravitational contribution is of the form 

Finally the curvature contribution to the potential energy 
is 

Thus the total potential energy 
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we get 



Each term in this expression is similar to the potential 
energy of a harmonic oscillator which on the average (in 
the classical limit) must equal k-QT/i so that 

A^{pg + -iq^^Kq^){\h{q)\')^kBT {A2) 

The power spectrum can be written asQ 
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Substituting (Al) in (A3) gives 
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Combining (A2) and (A4) gives 
C{q)^ 
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Since the potential energy is a quadratic function of 
the dynamical variable ft.(q), the distribution of heights 
Ph — {5{h — h{x.))) will be Gaussian. The difference 
Ah between the highest and lowest point on a randomly 
rough surface with a Gaussian height distribution, 

P{h) - (2^/.2„J-l/2cxp(-/lV2/^rms), 

depend on the size (area Aq) of the surface, but is in 
general much higher than the root-mean-square ampli- 
tude /irms- For example, if the height correlation function 
C(|x|) = {h{x.)h{0)) has a characteristic decay length 
C (i.e., C(O/C(0) « 1) then there will be roughly 
N = Aq/£^'^ uncorrelated points on the surface. Thus 
Ah w 2hi where hi is obtained from Jj^ dhP{h) ^ 1/N. 
For >> 1 this gives 
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For example, if A — 1 cm^ and ^ = 1/xm we get N = 10® 
and hi « 6/irms and thus Ah « 12ft,rms. If instead 
N — 10'' we get Ah « 8/irms, so that Ah is not very sen- 
sible to the exact value of N and in general Ah « lOft-rms- 
Restagno et al studied a 10 ^m x lO^m Pyrex glass sur- 
face region using AFM {N — 512 x 512 data points) and 
observed a hight fluctuation to be about 2 nm. This is 
about 10 times higher than the root-mean-square height 



fluctuation /irms expected for frozen capillary waves, and 
is roughly what is expected theoretically in this case. 

Appendix B: Stress distribution in adhesive 
contact between smooth spherical surfaces 

Greenwood and Johnson 30J have presented a very sim- 
ple and useful approach to the contact mechanics be- 
tween solids with locally spherical (or quadratical) sur- 
faces. They observed that the difference between two 
Hertzian pressure distributions, corresponding to two dif- 
ferent values a and b of the radius of the circular contact 
area, describes adhesive contact between the solids. The 
stress for r < a 



a = A{a^-r^)'^^-B{b^- 
where A > B > and b > a, while for a < r < 6 

^ =_i3 (52 _ ^2)1/2 

SO that (7 is continuous for r — a. In this model the 
spheres are in contact for r < a. For r > a no contact 
occur but there is an attractive interaction between the 
walls for a < r < b, while for r > 6 no forces act between 
the walls. Thus, we can consider the circular strip a < 
r < 6 as the crack tip process zone. Since the stress is 
finite everywhere on the solid surfaces, a nonzero crack 
process zone is necessary in order for the work of adhesion 
7 to be nonzero. If we introduce t] — b/a and denote the 
maximum compressional stress (which occur for r — 0) 
as (7max7 and the maximal tensal stress (which occur for 
r = a) with cja, and if we measure r in units of a we get 
for < r < 1 



.1/2 



(7/2 - 1)1/2 



(1-r^) 



2x1/2 



(772-1)1/2^'^ ^ 

while for 1 < r < 77 

.2 _ 2x1/2 



(7,2 



1)1/2^^ ) 



The applied normal force or load is given by 

2tt 



r<b 



and the normal force acting in the contact area 
Fo = / d^x a(x) = + ^^^(772 - 1) 

Jr<a •J 

Note that Fq > Fn as it must because the normal force 
in the contact area must equal the external load plus the 
additional "adhesion load" coming from the attraction 
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FIG. 31: The stress distribution P{ij) in the area of contact 
between two elastic spheres in adhesive contact and with- 
out adhesion (Hertz case). In the calculation we have as- 
sumed (Tmax = 1 GPa and the adhesive detachment stress 
ffa = —0.3 GPa, and the parameter 77 = 1.1. 



between the solids in the non-contacting (r > a) region. 
The pressure distribution in the contact area 



(ff'x 5{a — ct(x)) 



r<a 



This function is shown in Fig. 1311 for dmax — 1 GPa, 
da = 0.3 GPa and 77 = 1.1 corresponding to a crack tip 
process zone of width b — a = 0.1a. Note that P{(t) vanish 
continuously at cr = — (Ta- hi fact, it is easy to show that 
close to the detachment stress a — —aa, 



P{a) 



2{a 



where 



For (Ta = 0, CTi = CTmax 

(23). 



(772-1)1/2 
we obtain the Hertz solution, Eq. 
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